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ABSTRACT 


This report discusses the theory, design, implementation and testing of a 
personal computer-based Multi-Frequency Modulation (MFM) packet commu- 
nications system. Transmitter/receiver programs provide software drivers for 
D,A and A,D boards and perform symbol encoding, modulating, demodulating 
and decoding. The design and construction of a polarity coincidence correlator for 
receiver packet synchronization is presented. Experimental results show that the 
implemented MFM communication system conforms to theoretical analysis with 
acceptable bit error. Results also show that MFM can be uniquely adapted to a 


specific channel. 


THESIS DISCLAIMER 


The reader is cautioned that computer programs developed in this research 
may not have been exercised for all cases of interest. While every effort has been 
made, within the time available, to ensure that the programs are free of compu- 
tational and logic errors, they cannot be considered validated. Any application 


of these programs without additional verification is at the risk of the user. 
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I. INTRODUCTION 


Reliability and flexibility are two fundamental advantages of digital commu- 
nications. They are achieved through the speed and processing power of digital 
integrated circuits. Multi-Frequency Modulation (MFM) uniquely harnesses this 
power found in the modern personal computer to encode, modulate, demodulate 
and decode digitally formatted information, for sample rates currently up to 
100,000 samples per second. With computer-to-computer communication links 
implemented in a variety of mediums, such as wire, optical fiber, radio frequency, 
or acoustic, MFM readily adapts to a given medium and can emulate most ex- 
isting signal modulation formats. 

The focus of this thesis is the implementation of a MFM packet communi- 
cations system using industry standard personal computers (PC). Because of the 
limited processing speed of the PC, modulation and demodulation of the packet 
is performed in longer than real time. During the course of this project, trans- 
mitter and receiver software were developed to process data for establishing a link 
between PCs using D A and A’D data acquisition boards. A correlator was de- 
signed and constructed to synchronize the receiver to the transmitted packets us- 
ing a synchronization baud at the beginning of the packets. 

Analysis of the system phase response led to signal modifications which re- 
sulted in improved signal quality. Signal-to-noise performance was evaluated on 
the system for various packet constructions. The scope of the report encompasses 
the theory of MFM. a description of its specific development and implementa- 


tion, and a summary of the evaluation of the system’s performance. 


Il. THEORY OF MULTI-FREQUENCY MODULATION 


A. PACKET CONSTRUCTION 
The MFM signal set consists of “packets” of multiple tones which are ampli- 
tude and‘or phase modulated. These tones are present simultaneously during a 
subinterval of the packet known as a baud. Packets can be located arbitrarily in 
the frequency spectrum and time as seen in Figure 1. 
The follow definitions are used in MFM [Ref. 1: pp. 5-6]: 
T: Packet length in seconds 
AT: Baud length in seconds 
Ky 


o,: Symbol set. @, is the phase of the k” tone in the /* baud 


Baud length in number of samples 


L: Number of baud per packet 
At: Time between samples in seconds 
f, = 1/At: Sampling or clock frequency for D’A and A’D conversion in Hz 
Af=1/AT: Frequency spacing between tones 
K: Number of MFM tones 


To ensure the packet can be uniquely represented by the &, samples, the 
Sampling Theorem states that f,, the sampling frequency, must be greater then 
twice the highest frequency in the signal set [Ref. 2: pp. 46-56]. Since f, = k,Af 
and Af is fixed by AT, the highest tone is 4,/2 —1. Therefore, the signal set can 
consist of an arbitrary selection of tones from dc to f(,/2 — Af 

The analog representation for the MFM signal during the /” baud is 


k,|2 


x(t) = > Ay cos2nkAft+ by), (I- NAT < tS IAT. (1) 
=) 


The first baud begins at t= 0 and the last baud, L, ends at T = LAT . Sampling 
x(t) at t=nAt, where wm is discrete time, and substituting Ay — 1) A7 and 
At = AT/k,, produces a sampled output sequence of 


(MJ 5 
symbol Ik 


eg 


harmonic k MRS 
harmonic |X. 
FE QNY— 


k AT=Af'4 


-—~—_ 





time 


(eee eee 
Figure 1. NMIFM Signal Packet (after Ref. 1: p. 3.) 





k,|2 


x(n) = d A cox( 2 


a digital signal sequence of length k, samples. The k, point Discrete Fourier 


Transform (DFT) of x,(7) is 





ot +n) Ces Kes (2) 


k,|2 
X(k’) = a + kA yfeleud(k’ — k) + e /Pu5(k! — (k, —k))}, O<k' Sk, -1.(3) 
k=0 
From (3), it can be seen that the upper half of X,(k’) is the complex conjugate 
image of the lower half, which is consistent with the symmetry property of a real 


sequence [Ref. 2: p. 402]. 


B. SIGNAL GENERATION AND DEMODULATION 

In the past, generation of x(t) was attempted by constructing multiple phase 
lock loops, which were harmonically related. In December 1985, LCDR Deborah 
E. DeFrank, then at NPS, designed a coherent multifrequency synthesizer as a 
first step in producing x,t) [Ref. 3]. This method proved to be difficult to im- 
plement due to the phase jitter in the phase lock loops and changes required in 
hardware in order to change the harmonic frequencies used in the tone set. 

However, using the properties of DFT’s, a baud of x) can be generated by 
the host transmit computer by loading the first half of a complex valued array 
(0 to k,/2 —1) with the magnitude and phase of the tones to be included. To en- 
sure xm) is real, the upper half, (k,/2) to (k, — 1), is loaded with the complex 
conjugate of the values in the first half of the array at the image harmonics. An 
example of the DFT of a real x(m) having a period of k,=16 and 3 tones is shown 
in Figure 2. The Inverse Discrete Fourier Transform (IDFT) generates the real 
discrete sequence, x7), Which is clocked out of the computer thru a D A con- 
verter, at f, samples per second. A signal packet is generated by L repetitions of 
the above process. 

Demodulation of MFM is the inverse of the signal generation process. The 
analog signal, x(t), is sampled at f, samples per second and converted to digital 
format with an A’D converter. The sampled values are loaded into the real 
components of a k, point complex array, with the imaginary component set to 
zero. A DFT is computed of the array to obtain the complex frequency repre- 
sentation of the sampled input. Since the upper half of the DFT is redundant 


information, only the lower half is retained for further processing. 


C. PROPERTIES OF MFM 
Several important properties of MFM signals are presented in this section. 


(1) Orthogonality of signal during a baud. In continuous time 


i ne (1/2)AzAT k=i ; 
: Xy(t)x,{t) -5 rae (4) 


k Re(X(k)) -Im(X(k)) 
0 0 0 

I 0 0 

2 0 o- 
3 XR3 X13 
4 XR4 X14 
5 XRS X15 
6 0 0 

ql 0 0 

8 0 0 

9 0 0 

10 0 0 

11 XRS -X15 
ip XR4 -XI4 
13 XR3 -X13 
14 0 0 

15 0 0 


Figure 2. Data Structure (after Ref. 1: p. 9.) 


and in discrete time 
=k, -1 
a (1/2)A2k, kai 
0 


k #i i 


x,(a)x;(0) = : 


n=0 


The derivation of the orthogonality of harmonically related signals is well known 
[Ref. 4: pp. 152-154]. Two advantages of an orthogonal signal set are: 1)the 
noise, assuming additive white Gaussian noiseecAWGN), affects each transmitted 
tone independently, thus simplifying statistical computations, and 2)no tone in a 
baud interferes with any of the other tones. This is shown by representing each 


tone in (1), neglecting phase and magnitude, as 


5,(t) = cos(2nkAft) rect(t/AT) (6) 
where 
; tL 0 a <AT 
fect (au) — (7) 
0 otherwise 


The Fourier transform of 5,(z) is 
Sf) = ATsinc(f — KAJAT. (8) 


The spectrum of three adjacent tones, plotted in Figure 3, shows the peak of each 
tone coinciding with zero crossings of the spectrum of all other, thus minimizing 
their mutual interference. 

(2) Autocorrelation function (acf). Assuming x(n) is periodic, 


the circular acf of (2) is 


el 


rp) =) x(n)x(n@p), OS p<ky I (9) 


n=0 


where @ is a left circular shift. For a white bandlimited sequence with an even 


number of harmonics, the acf is given by 


sin(zKp/k,) 


r,(p) = (1/2)A7k, cos(2nkgp/k,) eS) 


(10) 
where ky = (k, + k,)/2, the midband harmonic, and K =k, —k,+ 1, the number 
of tones in the baud. Figure 4 shows the acf of a bandpass sequence with 
k, = 236, k, = 68 and &; = 83(K =16, 4, = 75.5). Mhese signal parametersecon- 
form to those of the synchronization baud, which will be discussed later. Note the 
peak of the acf occurs at p(0), this feature is fundamental to synchronizing the 
receiver to the incoming packet. 

(3) Matched filter. It is known that a matched filter maximizes 
the signal-to-noise ratio for additive white noise [Ref. 4: pp. 88-89]. The DFT 
of x(7) at the frequency &k is identical to the output of a filter matched to x,(z). 
This is due to the orthogonality property of the MFM signal and the linearity of 
the Dri [Retort 7-21. 


D. MODULATION 


Modulation is the process of encoding the source information onto a bandpass signal 
with a carrier frequency f.. This bandpass signal is called the modtlated signal s(1), 
and the baseband source signal is called the modulating signal m(1) (Ref. 5: p. 204]. 


L 


OAs “Cree 1.0 


2 


SP eG Vevsie MANU Uli 
0.0 


AIG 


68.0 69.0 70.0 720 72.0 7310 
HC INes. 


Figure 3. Orthogonal tone spacing. 


-0.4 -0.2 


The MFM signal can uniquely accomplish modulation in a number of ways. 
The signal x,(7), as defined in (2), can be modulated in amplitude, frequency, and 
phase, by translating the message into changes in A,,k and ¢, respectively. Note, 
any combination of the modulation types is also possible, such as changing am- 
plitude and phase to produce quadrature amplitude modulation (QAM). The 
signal, s(t), is called a bandpass signal. However, MFM can be bandpass or 
baseband and through multiplication with a carrier frequency x,(7) can be trans- 
lated to any frequency band desired. The signal sets considered in this thesis are 
bandpass and modulated using quadrature phase shift keying (QPSK) and dif- 
ferential quadrature phase shift keying (DQPSK). 


1870" 12.07 14,0 1620 


8.0 


MAG(R(P) } 


£.0 





a0 220 






a2 64.0 95.0 128.0 160.0 192.0 224.0 256.0 
' ie 
Figure 4. ACF of a white bandpass sequence (after Ref. 1: p. 16.) 


(1) QPSK. Conventional QPSK converts a digital input into 
four modulation voltage levels(symbols) to determine the phase of the transmitter 
output. To minimize the probability of symbol error the phases are spaced at 
multiples of 2/2 . A plot of the complex envelope of one tone of xm) is shown in 
Figure 5. The angle @ can take on the values of + 2/4 and + 37/4. 

Encoding MFM with QPSK is accomplished by loading the 
complex frequency domain array with the appropriate phase information. For 
example, if the digital input is ’0110...’. The first symbol, ‘01’, would be loaded 
into the frequency bin, k, as Re[.X(k)] = —A,//2 and Im[X(k)] = A,//2 , where 
A, is the amplitude of tone k. Likewise, the second symbol ‘10’ would be loaded 
into the second bin, k+1, as Re{X(kK+1)J= Vis D and 


Imaginary 
(quadrature) 


‘Ol 






Real 
(in phase) 


en 


Figure 5. QPSK signal constellation. 


Im[A(k + 1)J = Bee) WO This would continue until all the tones in the baud were 
filled or there were no more symbols. As mentioned earlier, the discrete time do- 
main signal is produced by taking the IDFT of the complex frequency array. 
Each successive baud is encoded similarly. 

Decoding QPSK directly into bits is accomplished easily as fol- 
lows. Assuming a coherent receiver, decoding requires evaluating the polarity of 
the real and imaginary components of each frequency bin. Notice in Figure 5 
that the symbol mapping uses Gray encoding. This reduces the probability of 
bit error because errors caused by AWGN are likely to occur when the adjacent 
symbol is selected for the transmitted symbol; thus, the symbol error will contain 
only one bit error. Gray encoding also allows decoding straight into bits, with the 
right bit determined by the polarity of the real component and the left bit by the 
polarity of the imaginary component. The digital signal is obtained by succes- 
sively decoding each tone of each baud. 

(2) DQPSK. QPSK requires strict phase coherent regeneration 
of the sampling frequency to eliminate phase ambiguity. This results in a complex 
Synchronization design Or distribution of the sampling clock frequency to each 


receiver. DQPSK provides a practical solution to the phase uncertainty problem 


Previous Symbol Previous Symbol 


‘00° “One 
00 Co 
01 10 
input 
symbol 
11 
ye next O01 V1 
symbol 
BOs 11 
Ol 
Jal, 10 11 
10 01 


00 00 
Figure 6. DQPSK encoding scheme. . 


but at a cost of approximately 2.3 dB more in SNR in order to obtain the same 
probability of bit error as coherent QPSK [Ref. 6: p. 442]. Encoding MFM with 
DQPSK is similar to QPSK; however, DQPSK translates the original symbol set 
into a second “differential” symbol set, which is then encoded as QPSK. Trans- 
lated symbols are generated based on the input symbol and the previous trans- 
lated symbol. Figure 6 shows this translation. Notice that, regardless of the 
previous symbol, an input of ‘00’ generates a new symbol in the same quadrant 
as the previous symbol. ’An input of ‘01’ rotates the new symbol + 2/2 radians 


from the previous, ‘!1’ rotates 2 radians, and “10’ rotates — 2/2. In the receiver, 


decoding is performed by determining the phase difference between successive 


pairs of tones, using complex arithmetic. The phase difference, Ad@, is found from 
gid? — eligi)" = oi poe eli ~JFi-1) (11) 
and in the absence of noise will be 0, + 2/2, or z radians. To realign the complex 


signal to the original constellations, e/4° is rotated by + 2/4 radians and the in- 
phase and quadrature components are decoded as with QPSK. 


Il. SYSTEM DEVELOPMENT 


Multi-Frequency Modulation joins the basic methods of digital communi- 
cations theory and signal processing. The theory behind MFM is not new [Ref. 
7], however, the ability to implement it inexpensively is linked to the recent de- 
velopment of inexpensive DSP chips for personal computers. The advent of 
packet switching in data communications also makes MFM a preferred choice for 
the MODEM because of its “packet” like format. 


A. SYSTEM DESCRIPTION 

An MFM communications system is shown in Figure 7. Information, de- 
noted by (mn), is modulated into a frequency band that will propagate over the 
available channel. The receiver converts the noise-corrupted signal, r(t), to an 
estimate of the source information, m(n). Theoretically, (m2) is any signal that 
can be represented in a digital format. To be practically useful, the conversion 
process must be real time. However, due to signal conversion speeds available for 
this research project, the input information had to be restricted to stored data 
files. The ideal transmitter output, x(f), is described for the /” baud by (1). 


n(t) 









information 









Transmitter Channel Receiver 


m(n) 


Figure 7. NIFM Communication System. 


B. BLOCK DESCRIPTION 
1. Transmitter 
The transmitter is subdivided into functional blocks, shown in Figure 8. 
As mentioned above, m(n) has been digitized and stored in a file. Therefore, the 


input is a string of binary digits. Processing is performed on a baud-by-baud 
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Figure 8. Transmitter Functional Block Diagram. 


basis, until the end of the data file, or maximum packet length has been obtained, 
whichever comes first. 

The encoder converts input symbols into complex values stored in the 
frequency domain array. The value of the symbols depends on the type of mod- 
ulation. For example, QPSK symbols are two bits long and are encoded as pre- 
viously discussed. The discrete signal, produced by computing the IDFT of the 
complex-valued frequency domain array, is loaded into the packet storage area. 
The controller determines the parameters of the packet based on the modulation 
type, baud, and message size; it then sequences the input data through the 
transmitter one baud at a time. Once the message has been processed, the entire 
Stored digital signal, x(), is transferred out at the selected rate, f, samples per 
second, through a D/A converter. Depending on the channel's frequency re- 
sponse, filtering of the output may be desired to remove the high frequency am- 
plitude discontinuities introduced by the D/A converter. 

2. Receiver 

The receiver, shown in Figure 9, demodulates the MFM signal by re- 
versing the transmitter process. Data acquisition is the process that samples the 
analog signal at f, samples per second and converts it to digital format with an 
A/D converter. Though not shown in Figure 9, filtering of the input is recom- 
mended to bandlimit input noise and to reject out-of-band interference. The 


converted data sequence is stored until all data is acquired. Again, this is due to 


time limitations in the signal processing algorithms available in the system used 
in this research task. The stored real values are accessed one baud at a time to 
perform a k, point DFT. The first half of the resulting complex values are de- 
coded to obtain the amplitude and phase modulation information. As in the 


transmitter, the controller sequences the data through the system. 
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out 


Figure 9. Receiver Functional Block Diagram. 


Notice the delay from input to output is the time it takes to acquire the 
entire packet and process the first baud. In contrast, real time processing must 
complete the k, point DFT in the same or less time than it takes to fill up one 
baud length buffer. Processing can alternate between buffers, and data will flow 
continuously through the system at the sampling rate. The storage requirements 
and processing time using f, = 61440 Hz, L =30, and k, =1024 samples are given 
in Table 1 for a system with a real time and non-real time DFT capability and 
for the non-real time system used in this research. The AT implementation 
processing time highlights the need for a high speed vector processor to perform 
the DFT computations. 

3. Synchronization 

In digital communications, various degrees of synchronization are re- 

quired. Typical degrees are: 
e Carrier synchronization, 
e Bit or symbol synchronization, 


e Word or frame synchronization. 


Table 1. STORAGE REQUIREMENTS AND PROCESSING TIME 


tion 
Storage required 2048 sans a, 
(# of samples) 2 


Time required for 
Slow > 16M 10 x 103 


processing 1024 
These have slightly different meanings and forms depending on the system. For 












samples (msec) 


the MFM system, carrier and frame synchronization may be required. 

Carrier synchronization is required if the modulation scheme is coherent, 
such as QPSK. This means generating an f, at the receiver in frequency and phase 
coherence with the transmitted f£. Though the packet does not contain the f, 
harmonic, every tone is harmonically related to it. The packet can be constructed 
to have a pilot tone separated from the modulation tones from which f, can be 
derived. For example, the modulation tones could range from k, to k,_, , giving a 
band of j + 1 consecutive tones. The pilot tone, k,, would then be placed several 
tones aWay to ensure minimum band interference from the message during its 
extraction. 

As mentioned previously, MFM, as implemented, is in the form of a 
transmission packet. To acquire the packet the receiver must know when to start 
sampling. This is accomplished by frame synchronization. Typically, unique 
words are inserted to mark the start of each frame. In MFM, the unique word is 
called a “synchronization baud”, and it is added at the beginning of each packet. 
This baud is generated similarly to other bauds, except the tones and phases are 
predetermined. Acquisition of the received signal starts after successful detection 
of the synchronization baud. [Ref. 5: pp. 511-512, 8: pp. 293-295] 


IV. SYSTEM IMPLEMENTATION 


In the previous chapters the tools and framework for developing MFM 


packet signals were established. One actual realization of this form of communi- 


cations will now be presented. Though the IBM Personal Computer (PC) was 


used, MFM can be implemented in a variety of ways. The PC, however, is ide- 


ally suited for MFM implementation for the following reasons: 


e Increased I/O channel maximum throughput rate of approximately 


100 KHz using direct memory access (DMA). 


Digital signal processing allows encoding, decoding, modulation, demodu- 
lation, and channel equalization. 


Signal processing algorithms are available in high level languages. 
Packet construction is easily modified to conform to various channels. 
External hardware is easily interfaced. 


Cost is low. 


Successful data acquisition and decoding of (1) was the first and most im- 


portant goal of this thesis. Other goals affecting software and hardware develop- 


ment were the following: 


A. 


Maximize packet size. 

Develop synchronization circuitry. 

Transmit. receive ASCII files. 

Develop software for testability and flexibility. 


Develop software for statistical testing. 


SIGNAL PARAMETERS 


Software and hardware implementations can be easily modified to receive any 


packet construction. However, all software conforms to the signal parameters in 


Appendix A, established between NPS and NOSC for an acoustic application. 


This is a bandpass signal in the band from 16-20 KHz using a clock rate of 


61440 Hz. All baud lengths are powers of two, allowing utilization of 


hardware software Fast Fourier Transform (FFT) algorithms to speed frequency 
transformation processing. 

Various baud sizes are available to give greater flexibility in adapting the 
packet to test specific channel parameters. In DQPSK, where the information is 
represented by the phase difference between adjacent tones, channel phase dis- 
tortion affects the shorter baud more, due to their larger Af’ The maximum 
tone-to-tone phase error introduced by a linear phase channel is shown in 
aie 2. 


Table 2. PHASE ERROR FROM A LINEAR CHANNEL 


Phase error 
(degrees) 





Obviously the longer baud would be preferred when differential coding is between 
adjacent tones. However, if the channel introduces time-related distortion, like 
noise bursts, electrical glitches, or propagation fluctuations, a short duration baud 
is desired. and differential coding should be between the same tones on adjacent 
bauds. 


B. SOFTWARE 
1. Transmitter 
Generation of the MFM signal has been accomplished in previous re- 
search at NPS. Utilizing previous hardware and core software, this project ex- 
panded the transmitter software to provide encoding of 16-QAM and QPSk, 
encoded data files using DQPSK, constructed maximum size packets, generated 
a synchronization baud, and provided greater flexibility in initializing the DMA. 


Actual transmitter hardware will not be discussed; it has been the subject of 
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previous research [Ref. 9]. Two new transmitter programs, TRANSMIT and 
XMITMES, will be discussed. 
aw TRANSMIT 

The program TRANSMIT provides maximum flexibility in con- 
structing and encoding signal packets. It is an excellent tool for basic system 
testing and can be used as a training aid to demonstrate various digital modu- 
lation schemes. A procedural flow diagram of TRANSMIT is shown in 
Figure 10. SelectBaud is a user interactive procedure for establishing the pa- 
rameters of a 16-20 KHz packet and for choosing the encoding scheme to be 
utilized. Parameters k&, and L are selected. From k,, the lower and upper 
bandlimits, k, and k,, are set. This generates a bandlimited packet with 
k, —k,+ 1 tones. The data storage requirement is k,L, representing the total 
number of samples in the packet which are clocked out. The encoding scheme 


determines the path for encoding the complex frequency domain array. 


ScaleBaud ComplexF FT SealeData DMAINIT 





EncodeQPSK 











EncodeQ AM 





Figure 10. TRANSMIT algorithm. 


EncodeQPSK begins by displaying a four symbol QPSK constellation 
used as a reference in selecting the symbols over the band. The symbols for the 
K tones in the baud can be selected in a variety ways: 

e Symbols for all tones are randomly selected from a random generator. 
¢ The symbol for each tone is selected by the user. 


e Individual tones may be removed from the band. 


This last feature allows construction of a baud with an arbitrary number of tones 
within the band determined by k, and k,. Selected symbols for the band are 
loaded into the complex frequency domain array with their complex conjugate 
image frequency. EncodeQAM is functionally the same as EncodeQPSK, except 
the symbols conform to a 16-QAM constellation [Ref. 1: p.27]. 

To obtain the real discrete time-domain sequence of the encoded 
baud, ComplexFFT computes the inverse FFT. ComplexFFT consumes the ma- 
jority of the processing time in the program due to the complex arithmetic oper- 
ations required, thus restricting the overall throughput of the system. For real 
time processing, the FFT algorithm must be accomplished by a hardware signal 
processor. 

Each value in the time domain sequence is represented as a real data 
type. occupying six bytes of memory. ScaleData converts these values down to a 
one byte format acceptable to the D/A converter and places them into a packet 
Storage buffer. EncodeQPSK, ComplexFFT, and ScaleData are executed for 
each baud, until all L baud have been processed. To transmit the packet out of 
the computer, DMAINIT transfers samples at f, samples per second over the 
DMA channel to D A converter [Ref. 9,10]. 

b, AMITMES 

To demonstrate the suitability of MFM for transferring information 
from a source to a sink, the program XMITMES was written to transmit an 
ASCII file encoded using DQPSK. Affixed to the beginning of the packet is a 
synchronization baud. As shown in Figure 11, XMITMES has an even simpler 
structure then TRANSMIT, because all tones in the band are encoded using 
DOQPSK. 

The synchronization baud is a predetermined sequence generated by 
SvncBaud. This baud is constructed as are all other baud, except k, is fixed at 
256, and tones 68 to 83 are encoded with the same random symbol pattern re- 
gardless of the packet construction or input message. This synchronization se- 
quence occupies the first 256 values in the packet buffer and therefore is the first 


to be clocked out of the computer. 
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Figure 11. XMMILTTNIES algorithm. 


As messages can be of various sizes, TailorPacket sets the maximum 
number of baud required for encoding. This is determined by dividing the num- 
ber of characters in the message by the number of characters that can be en- 
coded. 

DiffEncode encodes the message file into the complex frequency do- 
main array. It reads one character at a time; then breaks the eight bit character 
into four 2-bit symbols. The symbols are DQPSK encoded and stored in the 
frequency array. Once encoded, processing and signal output by ComplexFFT, 
ScaleData, and DMAINIT are the same as in TRANSMIT. 

2. Receiver 

Like the transmitter, the receiver requires hardware to interface the com- 
puter with the channel and processing software to demodulate the MFM signal 
packet. The channel interface is a high-speed data acquisition board, model 
DASH-16F, manufactured by MetraByte Corporation [Ref. 11]. Also at the re- 
ceiver, synchronization circuitry is required to detect the beginning of the packet. 
The synchronization circuitry is discussed in Section C. Implementation of the 
processing software assumes synchronization and a receiver f, the same as at the 


transmitter. Two receiver programs, RECEIVE and RECMES will be discussed. 
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a. “RECEIVE 

The program RECEIVE, shown in Figure 12, processes and displays 
only the first baud of a packet. It is an excellent tool to experiment with the 
DASH-I6F’s software drivers and it provides a quick indication of the overall 
system performance. Two procedures are available for data input, GetData and 
AcquireData. GetData reads in the time domain sequence file generated in the 
transmitter, eliminating all data acquisition hardware in the transmitter and re- 
ceiver; thus ensuring perfect synchronization and accurate evaluation of all en- 
coding, decoding, and processing software. The computational SNR of the system 
is determined under this condition, as the only noise present is from round off 


and truncation errors introduced in processing. 


Figure 12. RECEIVE algorithm. 











ComplexFFT 


GraphData 






Analog data acquisition is performed by the procedure Acquiredata. 
It initializes and controls the DASH-16F using procedures written by Quinn- 
Curtis [Ref. 12]. AcquireData allocates memory to store the sampled values 
transferred from the board using the DMA controller. Direct Memory Access is 
the only data transfer mode capable of transferring data to memory at the re- 
quired f,, without disruption by other interrupt processes in the computer. Other 
important initialization parameters are triggering source and the number of sam- 
ples to be collected. The A/D may be triggered from two sources, a programmable 
interval timer or an external trigger source. The programmable interval timer 


divides either a 1 MHz or 10 MHz clock to derive the sampling rate of the trig- 


21 


ger. Since this method cannot produce an arbitrary f,, external triggering is used. 
After initialization, conversions take place on the positive transition of every 
trigger until the specified number of samples have been acquired and transferred 
to memory. Collected data is then converted into a format acceptable for further 
processing. 

DiffDecode determines the encoded symbols by differentially decod- 
ing the complex frequency array transformed from the sampled data by 
ComplexFFT. The output of the program is the frequency spectrum plot in Fig- 
ure 13. The lower plot is the spectral response of the k,/2 tones in the baud; the 
upper plot represents only tones from k, to k, . The color on the original display 
indicates the phase quadrant of a given tone. Graphing the data in this fashion 
provides quick qualitative analysis of the frequencies and their phase informa- 


tion. 





Figure 13. Baud magnitude spectrum. 


b. RECMES 
RECMES, shown in Figure 14, demodulates the ASCII encoded 
transmission produced by XMITMES. It differs from RECEIVE in that it can 
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Figure 14. RECMES algorithm. 


process a multiple baud packet. The user interactive procedure, PacketSetUp, 
tailors processing to the expected receive packet, using the inputs k, and L. 
AcquireData samples and stores a memory segment of data regardless of the 
packet size, Lk, However, once stored, ConvertData, ComplexFFT, and 
DiffDecode process only L baud of the data. DiffDecode combines four 
differentially decoded symbols into one byte, representing the ordinal number of 
an ASCII character. To reconstruct the message, the characters are transferred 
to text file MESSAGE.DAT until processing is complete. For convenience, 


ShowMessage displays the recovered message. 


C. HARDWARE 

Synchronization of DQPSK MFM is obtained from a hardware correlator 
‘that is external to the host receiver computer. The 128 point correlator, illus- 
trated in Figure 15, provides the data acquisition board with sampling triggers 
synchronized with respect to time of arrival of the packet. Using only the polarity 
of the analog input, it functions as a matched filter to the last half of the 256 
point synchronization baud. This type of correlator is referred to as a polarity 
coincidence correlator(PCC). 

The hard limiter used to obtain the polarity information in the analog input 
is a fast, high precision, high gain, operational amplifier. During positive and 
negative portions of the input the output is +5 Vdc and O Vdc respectively. 
This unipolar signal is synchronized to the receiver's f, as it is clocked through a 
128 point serial shift register. To minimize bit instability, due to the shift register 


input being asynchronous to f,, the hard limiter slew rate should be as large as 
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Figure 15, Polarity coincidence correlator. \ 


possible. For example, during testing, 15 bits out of 128 were unstable using an 
MC1747 general purpose operational amplifier with a slew rate of 0.7 V/psec. 
This was reduced to 2 bits using the faster LM31I8 with a slew rate’ of 
70 Vipsec. 

The 1’s and 0’s buses store the time-reversed polarity sequence of the syn- 
chronization baud. The voltage on the buses represent the correlation between the 
stored sequence and the sequence in the shift register. As the sequences come into 
alignment, the 1’s bus voltage increases while the 0’s decreases. This inverse 
property is combined by a differential amplifier to give the total correlation. 
When both sequences match, the 1’s bus voltage is +5 Vdc and the 0’s is 0 Vdc, 
giving a maximum differential voltage of 5 Vdc into the threshold detector. The 
threshold detector generates a synchronization trigger on detection of the corre- 
lation peak. A latch is set, enabling the D flip-flop to pass sampling triggers for 
data acquisition at the packet’s beginning as required for demodulation. 

During design and testing the correlator progressed from a 16 to 128 point 
shift register. With each shift register output bus connected through a resistor 
voltage divider network, progressive testing was necessary to ensure power sup- 
ply current ratings were not exceeded. Maximal-length sequences Were used as 
test correlation sequences due to their unique two-valued autocorrelation func- 


tions that are easily determined by finite-field arithmetic [Ref. 8: pp. 368-375]. 
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V. PERFORMANCE EVALUATION 


Initial testing of the MFM system was concerned with quickly determining 
the quality of the demodulated signals. Figure 13 provided an excellent repre- 
sentation of the demodulated and decoded signal and was used throughout sys- 
tem development. This plot uniquely displays the frequency domain magnitude 
and phase. Since QPSK symbols are at least 90 degrees apart each color (on the 
original display) represents a specific quadrant in the complex frequency plot. For 
example, a green line indicates a tone whose phase is in the second quadrant and 
whose magnitude is represented by the length of the line. Although each color 
represents a broad decision region in phase, it is still clear a phase shift of at least 
180 degrees occurs across the 4 KHz band. This shift is primarily due to the linear 
phase Butterworth filter used to simulate the system channel, and it contributed 


to the decision to develop a DQPSK encoded signal. 


A. SYSTEM FREQUENCY RESPONSE 

After development of all hardware and software, testing revealed significant 
phase fluctuations across the frequency band. The system phase response, plotted 
in Figure 16, was determined by subtracting the transmitted phase from the re- 
ceived phase for each tone. Ideally the response should be smooth. With 
DQPSk encoded from tone-to-tone, phase differences introduced by the system 
will result in reduced tolerance to additive noise. For example, notice the phase 
difference of 0.4 radians between tones 296 and 297. When noise is added, de- 
coding the symbol in error will occur at a lower noise level than it would if the 
phase difference was zero; thus a higher SNR is required for low error rate de- 
coding when the channel introduces tone-to-tone phase fluctuations. The tone- 
to-tone phase fluctuations were substantially reduced, as shown in Figure 17, by 
decreasing the magnitude of the encoded tones in the transmitter which moved 
the operating range of the D/A converter to a more linear region. Further anal- 
ysis of the effect of the D/A and A/D converters on phase shift is recommended 


in order to be able to specify them properly in future designs. 
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Figure 16. Initial system response. 


A second important feature of the system phase response is its average slope 
across the band. The phase slope is affected by synchronization timing as shown 
in Figure 18. When the sequence through the correlator was delayed 2 bits, or 
2/61440 = 32.55 psecs, relative to the analog input to the receiver, thus delaying 
data acquisition accordingly, the phase slope changes as predicted by the time 


delay Fourier transform theorem; 
(Cee > Mien ae (12) 


This shows the phase slope in radians for a signal delayed by 7, seconds, is 2nT, 
radians/Hz. For ‘example, at f= Afk; = (60)(300) = 18000 Hz and 
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Figure 17. Improved system response. 


T, = 32.55 psecs, the phase shift is 3.7 radians. This compares closely to the ac- 
tual phase shift of tone 300 in Figure 17 and Figure 18. 

In summary, system frequency response analysis identified phase shifts in- 
troduced by the system hardware and by the synchronization timing. The mag- 
nitude of the encoded tones and the synchronization sequence delay were selected | 
experimentally in order to minimize tone-to-tone phase fluctuations and provide 


the comparatively flat response shown in Figure 18 for all subsequent testing. 
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Figure 18. System response with 2-bit synchronization delay. 


B. SNR PERFORMANCE 

The performance of the system when corrupted by additive noise is deter- 
mined by its output signal-to-noise ratio (SNR). For QPSK-encoded MFM the 
mean of the real and imaginary parts of the 2K coefficients of the DFT represent 
the reccived signal amplitude, and their variance represents the noise power. The 
output SNR is defined as the ratio of the square of the mean to the variance of 
each for these 2K coefficients [Ref. Ll: p. 25]. 

Figure 19 shows the system SNR, which is the output SNR when there is no 
additive noise in the channel, versus the frequency spacing, Af, of the five baud 
types. As would be expected, DQPSK performs better with a smaller Af because 


the phase difference between adjacent tones is smaller when Af is smaller. 
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Figure 19. Nilaximum SNR output. 


In Figure 20, the SNR out of the receiver is shown versus input SNR which 
was varied by adding noise to the analog signal. Theoretically, the output SNR 
should equal the narrowband input SNR in QPSK [Ref. |: pp. 24-25]. The re- 
duction in output SNR at higher input SNRs is due to each baud approaching 
its maximum output SNR as set by system noise and shown in Figure 19. 
Table 3 lists the bit errors at various SNR input levels. There were approxi- 
mately 2500 bits transmitted for each entry. This clearly indicates the worsening 
effect of the channel on bauds with greater tone spacing Af at high SNRs. 

A direct comparison of the MFM signal with commercial modems is difficult 
due to the wide variety ‘of techniques that are used to achieve high-speed data 


transfer. However, with a bandwidth efficiency of 2 bits/s/Hz, a DOQPSK MFM 
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Figure 20. SNR performance. 


signal can achieve a throughput of 6000 bits per second on a standard switched 
telephone line having a bandwidth of 3 KHz. Though not implemented in this 
project, several techniques could be utilized to increase the bandwidth efficiency, 
such as increasing the constellation size and/or using data compression algorithms 
as in commercial modems [Ref. 13, 14]. Using these techniques MFM could 
achieve a bandwidth efficiency of 4-8 Bits/s/Hz, which compares favorably to 


current high speed modems. 
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Table 3. BIT ERRORS IN 2500 BITS TRANSMITTED VS BAUD TYPE AND 


SNR. 
SNR (dB) 
ee ee eee eee 
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VI. CONCLUSIONS AND FURTHER STUDY 


MFM is well suited to the signal processing environment of the personal 
computer. As implemented in this thesis, the packet signal provides a high speed, 
low bit error data transfer link between two industry standard computers. Syn- 
chronization hardware allows the link to use asynchronous data. 

The theory and properties of MFM have been discussed for a white bandpass 
packet signal. Packet construction and various modulation formats have been 
shown to be easily adapted to a given channel. The system was developed and 
implemented for testability and to be easily modified to accommodate a variety 
of applications. However, the data storage requirements and signal conversion 
speeds were unique to this project due to available hardware. As implemented, 
this system would require a multiple packet transmission to achieve practical op- 
eration. 

Analysis of the system’s phase response led to adjustments in the encoded 
tone magnitudes and synchronization timing. This substantially improved the 
svstem’s phase response. SNR test results indicate a smaller Af baud has superior 
performance in a linear phase channel. For Af less than 60 Hz, bit errors are 
acceptably low for input SNRs greater than approximately 15 dB. Exhaustive 
testing is required to compare bit error performance with conventional digital 
modulation schemes. 

Areas of further study should focus toward increasing the system bandwidth 
efficiency and implementation of hardware signal processing for real-time opera- 
tion. In fact, current research at NPS is directed at increasing the throughput 
speed of the system described. Improved bandwidth efficiency can be obtained 
through data compression algorithms and developing receiver software to decode 
the 16-QAM encoded signal generated by TRANSMIT. Finally to provide a 
greater flexibility in adapting MFM to channels with significant uncompensated 
phase distortion, the option should be developed to encode DQPSK on a baud- 
to-baud basis. 


APPENDIX A. DESIGN PARAMETERS 


Table 4. DESIGN PARAMETERS FOR A 1/15TH SECOND SIGNAL PACKET 
IN A 16-20KHZ BANDPASS CHANNEL. 
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APPENDIX B. OPERATING INSTRUCTIONS 


A. PRELIMINARIES 
1. Hardware setup 
a. System setup 
Setup system as shown in Figure 21]. All signal/trigger connections should 


be accompanied with a ground lead to minimize interference. 
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Figure 2t. System interconnection diagram. 


b. DASH-16F switch setting: 
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GAIN (12.5y) ----------2-------22-2- 22 nnee cee dipswitch A, B--ON; C, D, USER--OFF 
Address (300h) -------------------------------- dipswitch 9, 8--OFF; 7, 6, 5, 4--ON 


c. Trigger Requirements 
The D/A converter board uses negative logic as shown in Figure 22. The 
pulse width, PW, must be less then one psec to ensure proper signal output. This is due 
to DMA request and acknowledge hardware. All hardware on the correlator and 


DASH-I6F is positive edge-triggered. The trigger magnitude should be approximately 
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Figure 22. Trigger specifications. 
d. Power Supplies 
+8Vde supply ---------------------------------- HP6216A 30v 500ma (or suitable substitute) 
—8Vde supply ---------------------------------- HP6216A 30v 500ma (or suitable substitute) 
+5Vdce supply ---------------------------------- HP6216A 30v 500ma (or suitable substitute) 


2. Software setup 


a. Transmitter 


1. Convert DMAINIT.ASM and DMASTOP.ASM_ to BINary _ files 
[Ref. 15: pp. 91-93}., 


2. Place files DMAINIT.BIN, DMASTOP.BIN, COMPFFT.INC, — and 
FF1I87B2.INC in the same directory as TRANSMIT.PAS AND XMITMES.PAS. 
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3. Compile to disk TRANSMIT.PAS and XMITMES.PAS. 


b. Receiver 
1. Compile TP4D16.PAS to disk to create TP4D16.TPU. 
2. Place TP4D16.TPU in the same directory as RECEIVE.PAS and RECMES.PAS. 
3. Compile to disk RECEIVE.PAS and RECMES.PAS. 


B. MESSAGE TRANSMISSION PROCEDURE 

. As necessary conduct preliminary equipment setup. 

. Reset correlator (removes triggers to DASH-16F). 

. Run XMITMES until “Ready to Transmit”. 

. Run RECMES (match baud size used in XMITMES). 


5. When receiver prompts “Ready to acquires”, press enter key on transmitter. 


mR WY N= 


6. Demodulated message is displayed on receiver screen and stored in file 
MESSAGE.DAT. 


TRANSMIT and RECEIVE are similarly executed. However TRANSMIT would 
require manual loading of the synchronization baud in Table 5 by individual loading 
each tone with the appropriate symbol. 


Table 5. SYNCHRONIZATION BAUD SYMBOL SEQUENCE 


| tone _| 68 | 69 | 70] 71 | 72 | 73 | 74 | 75 | 76 | 77 | 78 | 79 | 80 | 81 | 82 | 83 





[symbol] 3 {3 {2} 2;atrt2qtit2}3f4aistit4]3 fil 
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APPENDIX C. TESTING PROCEDURES 


A. RESPONSE TESTING PROCEDURE 
1. Conduct preliminary equipment setup. 
2. Run XMITMES and exit. 


3. Transfer file XMITDAT.DAT to the receiver computer. (Not required for subse- 
quent response testing if transmitted message and packet construction are un- 
changed.) 


4. Transmit message IAW message transmission procedure. (RECMES generates 
output file RECSTAT.DAT used by RESPONSE and STATS.) 


5. Run RESPONSE. (Ensure baud size is correctly edited into RESPONSE.) 
6. Graph output file RESPbaud.DAT. 


B. SNR TESTING PROCEDURE 


1. Conduct preliminary equipment setup. 


2. Add random noise source and true RMS voltmeter as shown in Figure 23. 
2.2KN 






Receiver 










True 
RNS 


Voltmeter 
Random 


Noise 
Generation 


Figure 23. Test equipment interconnection. 
3. Set SNR level. 


4. Conduct steps 1-3 in the Response testing procedure. 
5. Run Statistics (results are in file STATbaud.DAT). 
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APPENDIX E. TRANSMIT 


program TRANSMIT; 

(*Transmits variable length packets up to 61440. The length is 
determined by BCSTARRAY and can be take to a full segment of memory. 
No syne baud is used. All external files are assumed to be in the same 
directory”) 


type 
TNvector = array[0..4095] of real; 
TNvectorPtr = aTNvector; 
BCSTARRAY = array[-28673..32767] of byte; (*stores output samples*) 


var 
kx, (*baud size*) 
kl, (*lowest tone in band*) 
ke, (*highest tone in band*) 
NUMBAUDS , (*number of bauds in packet*) 
BAUDCOUNT, (*current baud being processed*) 
BYTECOUNT, (*number of bytes to be transferred by DMA*) 
I : integer; 
XREAL, (*complex freq/time domain arrays*) 
XIMAG > TNvectorPtr; 
QAM, (*encoding scheme) 
INVERSE : boolean; (“direction of FFT*) 
ERROR : byte; (*status of ComplexFFT*) 
BCST : BCSTARRAY; (*packet storage buffer*) 
THEFILE :file of byte; (“time sequence output file*) 
ANSWER : char; (*input variable”) 


($1 FFT87B2. INC*) 
(*SI COMPFFT. INC*) 


procedure SelectBaud; 
(*SelectBaud establishes QPSK or QAM encoding kx, kl, k2, and the number 
of baud in the packet*) 


var 
ANSWER > integer; 
MODANSWER- : char; 


begin 
window(20,10,80,20); 


(*Select encoding scheme*) 


repeat 
clirsienr; 
writeln('Select modulation desired'); 
writeln; 
writeln(' A --- QPSK'); 
writeln(' B --- 16 QAM'); 
readln(MODANSWER); 
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until MODANSWER in ['a oy Bs 
if MODANSWER in ['b','B ons 


QAM := true; 


(*Select buad size*) 
kx: =0; 
repeat 
clrscr; 
if kx < 0 then writeln('TRY AGAIN"); 
writeln('What is the length of the bauds (kx)?'); 
Wrieeln( i.e, 256, 512, 1024, 2048, 4096'); 
readln( ANSWER); 
case ANSWER of 
ANOR \koee SSDS 
512: kx:=512; 
1024: kx:=1024; 
2048: kx: =2048; 
4096: kx: =4096; 
end; (*case ANSWER*) 
if kx = 0 then kx := -1; 
mele akexee =O} 


(*Set bandlimits*) 
case kx of 
256: begin 
k1:=68; k2:=83; 
end; 
512: begin 
k1:=135; k2:=166; 
end; 
1024: begin 
alk — 26 Smke2 oo 
end; 
2048: begin 
k1:=537; k2:=664; 
end; 
4096: begin 
k1:=1073; k2:=1328; 
end; 
end; (*case kx*) 


(*Select number of baud*) 
write('How many bauds do your desire to transmit? '); 
readln(NUMBAUDS); 
window(1,1,80,25); 

end; (*SelectBaud*) 


procedure DisplayQAM; 

(*DisplayQAM shows the 16-QAM constellation in the upper portion of the 
screen. A window is set at the bottom of the screen for further 
interaction with the user. *) 


begin 


@ lliaseses 
SOLOsY« il, 1). 


4] 


write(chr(218)); 
for I:= 12 to 64 do 
write(chr(196)); 
(Wesel elasel( iL) a1) 3) 
for I:= 2 to 19 do 
begin 
GotoXY(11,1); write(chr(179)); 
GotoXY(65,1); write(chr(179)); 
end; 
GotoXY(11, 20); 
write(chr(192)); 
for I := 12 to 64 do 
write(chr(196)); 
write(chr(217)); 


window(12,2,63,19); 

write('This program encodes a 16 QAM multifrequency'); 
Weitein( sienale Je 

writeln('The vector space is shown for one frequency. '); 


writeln(' ,  chEGi79 ye 
writeln(' * * "  chs@iggyn. * * Se 
writeln(' 4 3 ',chr(179),! iz il se 
writeln(' 'chr(179),' ')3 
writeln( * * ,chr(179), | * * 23 
writeln( 8 7 5 lone Il 7/S))) . 6 5 3 
writeln( chicciy > ee ); 
Be In ; ; pees 2, ‘ 
writeln 7 7 Glin : * 7 
writein(' 12 pl ene , 10 9 »: 
La ay ; Dei , , pe 
writeln * fe ecelniae ; 3 r) : 
writein(' 16 iS " chrCig oye. 14 As ae 
writeln(' " cnr(i79), ye 

(*active window*) 
WeliclOw (elipez lis OneZon): 

end; (*DisplayQAM*) 

Ce SOROS Se oO Se oS eC o SSS sooo oSe Soc oeS Ss Seco o oS SSeS SSeS Se scooeSeSoSe soee soo seecece *) 

(* eee we) 


procedure DisplayQPSk; 

(*DisplayQPSK shows the QPSK constellation in the upper portion of the 
screen. A window is set at the bottom of the screen for further 
interaction with the user. *) 


begin 
@IGewiek 
eqoeepyOd ial g ie 
(re aeel( elnvet( 2 iis) je 
for I:= 12 to 64 do 
write(chr(196)); 
peaeat( elnset( 1S) )) je 
cepe Ie 2 ite ile) oo 
begin 
GotoXY(11,1); write(chr(179)); 
GotoxY(65,1); write(Cehr(179)); 
end; 


GotoXY( 11,20); 

write(chr(192)); 

ieee IC BES 1 sere) (5x5 Tolle 
write(chr(196)); 

write(chr(217)); 


wimdow( 12).2,63,,19)- 
writeln('This program encodes a QPSK multifrequency signal.'); 
writeln('The phase are shown for one frequency. '); 


writeln(' yeche( eran OE 
writeln(' * Nehr( 179), | ve oN 
writeln(' 2 ",chr(179),' il ey: 
writeln(' echrC Zone oF 
writeln( ,chr(179), ys 
writeln( ,chr( 179), ae 
writeln( | »chr(179), 3 
writeln( ,chr(179), oe 
writeln( | ,echr(179),_ 33 
writeln( ,ehr( 179), 3 
writeln( i @haset( IL 7/8) )S . ye 
writeln(' * Penn 1 7o0n # “OE 
writeln(' | Schr Gl7ojee 4 ye 
writeln(' Echroi79 ye ye 

(*active window*) 
wangowel2,21,80),25); 

end; (*DisplayQPSkK*) 

(* ee ee ee a a ee eee *) 

Wee ee mem ww mmm www www ew ewe ew ew we ee we eee ee meee eww eee em ween ewe wm eee ee coce *) 


procedure SelectQAM; 
(*SelectQAM selects and encodes the symbols for each tone in the band set 
by kl and k2. All symbols can be randomly or manually selected. *) 


var 
RESPONSE : char; 
ANSWER, I : integer; 
VECTORARRAY : array[0..4095] of integer; 

( Ww wm ww ww ww ww ow = www ow ww ww ww ow ow ww wn ww ww ww ww ww ow we ee we ee ee ee ee eee *) 


procedure EncodeData; 
(*EncodeData loads the frequency domain arrays XREAL and XIMAG with the 
QAM symbols in VECTORARRAY and their complex conjugate image. *) 


var 
J : integer; 
Wa Wael B saya lls 


begin 
fillchar(XREAL« ,sizeof( XREALA~) ,0); 
fillchar( XIMAG., ,sizeof( XIMAG~) ,0); 
igoie Jee Jeil tee Ike ele 


begin 
case VECTORARRAY[ J] o 
0 : TEMPI:= 0.0; 
Ion Snag ss Sie). Os 
Soo 9 Wawletec Sel, Ce 
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9..12 : TEMPI: =-30. 0; 
13.267 PEMPT: =-907 0; 
end; (*case VECTORARRAY*) 

case VECTORARRAY[J] of 


0 : TEMPR: = 0. 0; 
1,5,9,13 : TEMPR:= 90.0; 
2,6,10,14: TEMPR: = 30.0; 


See lS; TEMPR:=-s00 0: 
4,8,12,16: TEMPR: =-90. 0; 
end; (*case VECTORARRAY*) 


(*Load complex conjugate image*) 
XREALa[ J]: = TEMPR; 
XREALA[ kx-J] : = TEMPR; 
XIMAGa[ J]: = TEMPT; 
XIMAGa[ kx-J]:= -TEMPI; 
end; (*J; =k1*) 
end; (*EncodeData*) 


writeln(' Select one of the following for baud ',BAUDCOUNT,'.'); 
writeln; 
writeln('R Randomly select all ',k2 - k1 + 1,' vectors’); 
writeln('I Individually select vectors'); 
read1n( RESPONSE); 
UneieRESPONSH inelaer . nee 2 | 
fillchar( VECTORARRAY, sizeof(VECTORARRAY), 0}; 
case RESPONSE of 
eee R eetor 1i— hietommemdc 
VECTORARRAY| I] : =random(16)+1; 
ive = begin 


I:= kl; 
while I <= k2 do 
begin 
repeat 

clrscr; 
writeln(' Vector Selection Menu'); 
writeln('# vector desired"); 
wrttein( | 17 random vector'); 
writeln('18 no vector’); 
write('19 no more tones fene  ,1, 0. 


readin (ANSWER); 
until ANSWER in [1..19]; 
case ANSWER of 
1. . 16: VECTORARRAY{ I] : 
17: VECTORARRAY[ I] : 
18:3; (*nop*) 


ANSWER; 
random( 16)+1; 


Oe he? 
end; (*case ANSWER*) 
I: =I+1; 


end; (*while I*) 
end; (*'i' Ree) 
end; (*case RESPONSE*) 
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EncodeData; 
end; (*SelectQAM*) 


procedure SelectQPSK; 
(*SelectQPSK selects and encodes the symbols for each tone in the band 
set by kl and k2. All symbols can be randomly or manually selected. *) 


var 
RESPONSE Mehialis: 
ANSWER, I : integer; 
PHASEARRAY : array[0..4095] of integer; 


procedure EncodeData; 
(*EncodeData loads the frequency domain arrays XREAL and XIMAG with the 
QPSK symbols in VECTORARRAY and their complex conjugate image. *) 


var 
di : integer; 
TEMPR,TEMPI : real; 


begin 
fillchar( XREAL« ,sizeof( XREAL») ,0); 
fillchar(XIMAG~ ,sizeof(XIMAG~) ,0); 
POC OmK2 do 
begin 
case PHASEARRAY|[ J] of 
0 SMe sloits=) 0), Oe 
2 vTENET: = Sono; 
3} fh TEMP] =-8 00: 
end; (**case PHASEARRAY*) 
case PHASEARRAY|[ J] of 
0 SEMPRA OPO) 
1,4 : TEMPR: = 80. 0; 
2S 2 Mae SiO), (02 
end; (*case PHASEARRAY*) 
XREAL«[ J]: = TEMPR; 
XREAL«[ kx-J]:= TEMPR; 
XIMAG,~[ J]: = TEMPI; 
IMAG A K<-J)2= =-TEMPI; 
end; (*J: =k1*) 
end; (*EncodeData") 


Ciltrsicn 
writeln(' Select one of the following for baud ',BAUDCOUNT,'.'); 
writeln; 
writeln('R Randomly select all ',k2 - kl + 1,' phases’); 
writeln( ‘I Individually select phases’ 3 
readln( RESPONSE); 
Une RESPONSE in ['r°,'R', a at ds 
fillchar(PHASEARRAY, 5 izeof( PHASEARRAY) , ODE 
case RESPONSE of 
foe kh “fore — kl to k2ede 


PHASEARRAY[ I] : =random(4)+1; 


ere 1c ‘begin 
I:= kl; 
while I <= k2 do 
begin 
repeat 

clrscr; 
writeln(' Phase Selection Menu’); 
writeln('# phase desired'); 
writeln('5 random phase' ); 
writeln('6 no phase'); 
write('7 no more tones Tone “.l a as 


readin (ANSWER); 
until ANSWER in [1..7]; 
case ANSWER of 


1. .4: PHASEARRAY[ I]: = ANSWER; 
5: PHASEARRAY[ I]: = random(4)+1; 
6: ; (*nop*) 
7: I:= k2; 
end; (*case ANSWER*) 
I:=I+1; 
end; (“while I*) 
end; (*'i','I'*) 
end; (*case RESPONSE*) 
EncodeData; 
end; (*SelectQPSK*) 
(* SS SOS Sp Se SOS OOO SOS DES Sp OO Ope mmo mS ooo Spo so oS Sess seep oe Sse aeaasoo5oen *) 


procedure ScaleData; 

(*ScaleData converts each real value in array XREAL down to a byte and 
stores the byte in packet storage buffer BCST. INDEX sets the location 
in the buffer of each byte*) 


var 
UNDER Jee Mir : integer; 
DATA : byte; 
begin 
for J := 0 to kx-l1 do 


begin 
TEMP := round(XREALa[ J] + 126); 
if TEMP < 0 then 
TEMP := 0; 
DATA := TEMP; 
INDEX := J + (BAUDCOUNT - 1) * kx - 28673; 
BCST[ INDEX] := DATA; 
write(THEFILE ,DATA); 
end; (*for J*) 
end; (*ScaleData*) 
( Tew www ww ww ww we a a a a a a a a a ew eee ee eee ww ee ee ee ee ee ee ee eee ee eee *) 
(We nee enn enn ee ee ee ne ee ee eee eee *) 
procedure Dmainit(var BCST : BCSTARRAY; 
BYTECOUNT : integer); 
(*Assembly language procedure used to initialize and unmask the DMA for 
data transfer. The source must be converted to a BIN file.*) 
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external 'DMAINITS. BIN'; 


procedure Dmastop; (*Masks DMA, stopping data transfer“) 
external 'DMASTOP. BIN'; 


( BS SOS OO SO I te Cea ee Pager ea a *) 
procedure TINPUT; 

(“This procedure reads in binary file, HOUTFAl1.DAT, into an array, 
performs an FFT using ComplexFFT and displays the results 

with GraphData*) 


var 
NUMPTS : integer; 

(* COD CODD e Seo Se Sooo Deseo o pO Do OCDOCSoCe Oop eo oe oS OOD O UES Oo OMOo Sooo peorsooooS *) 

procedure GetData; 

(*GetData reads time domain sequence from file HOUTFA1. DAT into XREAL*) 


var 
WHERILE :file of byte: 
DATA > byte; 
dis 
GENE : integer; 
rdata : real; 

begin 
assign(THEFILE, 'HOUTFA1. DAT’ ); 
reset(THEFILE); 
new( XREAL); 
new(XIMAG); 


fillchar(XREALA,sizeof (XREALj),0); 

fillchar(XIMAGa ,sizeof (XIMAGa),0); 

Ota a Om OmKx—slNido 

begin 
read(THEFILE,DATA);(*read file one byte at a time*) 
TEMP := DATA; 
temp := temp and 255; (*reduces number of bits used*) 
(*to represent the input data”) 


KREALA[ I] := (TEMP - 126); (*loads XREAL array*) 
end; (*for I*) 
close(THEFILE); 
end; (*GetData™) 


Same we eee eee wwe Bee ewe eee Bee ee ee wee eww eB ew we Bee eB eee eee wee tee see ee ee He vs ) 


procedure PutData; 
(“Writes frequency domain arrays XREAL and XIMAG to file TFFTOUT. OUT*) 


var 
I : integer; 
OUTFILE : text; 


begin 
assign(OUTFILE, 'TFFTOUT. DAT"); 
rewrite(OUTFILE); 
for I := 0 to kx-1 do 
Writeln(OUIFILE,1,‘ ',XREALA)@)I,’ ‘*,XIMAGAf 1} ); 
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close(OUTFILE); 
end; (*PutData*) 


(* wi a a a a a a a a a a no i we a ee ee ee ee eee ee eee *) 


procedure GraphData; 


(*GraphData plots the kx/2 points of complex array XREAL and XIMAG. 
Symbols are decoded and represented by the color of each tones 


spectral line.*) 


var 
bee 
anyon 
COLOR, 
evan 
Y,YZ 


(*loop indices*) 


: integer; 


begin 
GraphColorMode; 
palette(2); 
Graphbackground(0); 
Gasca w (5 One nO nr Onn leOrele) i 
draw(30,170,300,170,1); 
clean 30, 155 300,155, je 
draw(30,80,300,80,1); 


(*every other dot equals one bin at 


draw(166,170,166,174,1); (x-axis 
draw(196,170,196,174,1); (*x-axis 
Giga @sOnelelone cmelelornle) |: (*y-axis 
ober 30), 25, 265255 UR (*y-axis 


IRS 2SR 


(*y-axis*) 
(*x-axis for full spectrum*) 


scale 
scale 
scale 
scale 


(*x-axis for zoom spectrum’) 


256%) 
marks* ) 
marks*) 
marks* ) 
marks*) 


(*Grouping of 4 horizontal lines to identify background color 


used for one of the symbols*) 

repeat 
GacawiG@oileelens OOmplin On)i: 
eleemn( Sil. Isril, SHO), aril , je 
draw(sit+2, 300,172, 2)5 
lieeuy( Sil, Iar8\,, SiO), Wars 32 

2 Slsriles 
BiaienIl WS 75g 
KXDIV2 := round(kx / 2); 


(*Decode symbols and assign color*) 
for I := 0 TO KXDIV2 do 
begin 


if (XREAL»[I] >= 0) and (XIMAGs[I] > 0) then 


COLOR := 0; 


if (XREAL«[I] < 0) and (XIMAGa[I] > 0) then 


COLOR := 1 


if (XREALa[I] < 0) and (XIMAGs[I] <= 0) then 


COLOR := 2; 





if (XREAL»[I] >= 0) and (XIMAGa[I] <= 0) then 


COLOR := 3; 


(*zoom spectrum*) 


if (I >= kl) and (I <= k2) then 


begin 

XZ := round(((I * 4096.0) / kx) - 1040); 

YZ := 80 - round(0. 5*(sqrt(sqr(XREALs[I]) + sqr(XIMAGa[I] )))); 
draw(XZ,80,XZ,YZ,COLOR); 

end; 


(*full spectrum*) 


X := round(I/KXDIV2 * 256 +30); 
Y :=170-round(20*1n(sqrt(sqr( XREALa[ I] +1)+ 
sqr(XIMAGa[ I] +1)))/1n(10)); 
draw(X,170,X,Y,COLOR); 
end; 
end; (*GraphData*) 


begin 
INVERSE := 
ERROR := 0; 
GetData; 
NUMPTS := kx; 
writeln('Performing FFT’); 
ComplexFFT (NUMPTS, INVERSE ,XREAL, XIMAG, ERROR); 
writeln('Error = ',ERROR,' hit the enter key’); readln; 


false; (*sets forward FFT*) 


PutData; 


GraphData; 
end; (*Tinput*) 
Cc BECO OD OOO ETE TN eC es ea ee ee *) 


begin 
repeat 
clrscr; 
assign(THEFILE, 'HOUTFA1. DAT' ); 
rewrite(THEFILE); 
new( XREAL); 
new( XIMAG); 
INVERSE: =true; 
QAM: =false; 


SelectBaud; 


if QAM then 
DisplayQAM 

else 
DisplayQPSkK; 


(*Packet contruction loop*) 
for baudcount := 1 to numbauds do 
begin 

if QAM then 
SelectQAM 

else 
SelectQPSK; 

writeln('Performing IFFT’); 
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ComplexFFT( kx, INVERSE , XREAL, XIMAG, ERROR); 
ScaleData; 
end; (*for BAUDCOUNT*) 


dispose(XREAL); 

dispose( XIMAG); 

close( THEFILE); 

BYTECOUNT := NUMBAUDS*kx-1; 

writeln('Press return to transmit'); readln; 


Dmainit(BCST,BYTECOUNT); 


repeat 
writeln('Transmit some more? (*yes or no*) "'); 
readln( ANSWER); 

Until IeANSWER in [| ns Neeeeves  ]: 

dmastop; 


window(1,1,80,25); 
until ANSWER in ['n','N']; 
Tinput; 
end. 
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APPENDIX F.) XMITMES 


program XMITMES; 

(*Transmits a syncbaud and message from file 'MESSAGE.DAT'. The message 
is encoded using QPSK. 'MESSAGE.DAT' is a text file. It should already 
exist before using this program. Output is used to collect data for 
TESTING*) 


const 
FIRST_ELEMENT = -28929; 


type 
TNvector = array[0..4095] of real; 
TNvectorPtr = aTNvector; 
(*Syne + 61440*) 
BCSTARRAY = array[ FIRST_ELEMENT. .32767] of byte; 


var 
kx 
Iealikee el Wie 
NUMBAUDS , MAXNUMBAUDS , 
BAUDCOUNT , BYTECOUNT, 
SYMBOLCOUNT , MAXNUMCHAR, 


MESSAGESIZE : integer; 
MAGNITUDE, 

CHARACTERS_PER_BAUD : real; 

XREAL, XIMAG :TNvectorPtr; 
INVERSE : boolean; 
TEMPBYTE , ERROR : byte; 

BCST : BCSTARRAY; 
eb 818 Tle :file of byte; 
Ways 1Ne IOs > text; 
ANSWER, 

YEXTCHAR : char; 


(=SlePrie7e2. INC } 
(=$1 COMPFFT. INC*) 


procedure SyncBaud; 
(*Process the synchronization baud and stores the 256 point time domain 
sequence at the beginning of the packet storage area. *) 


var 
SYNCREAL, SYNCIMAG : TNvectorPtr; 
Vy Weste : integer; 
SYNCDATA i byte; 
SYNCMAG ; real; 

begin 
new(SYNCREAL); 
new(SYNCIMAG); 


SYNCMAG: = MAGNITUDE; 
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(*load synchronization symbols*) 
fillchar( SYNCREAL~ ,sizeof(SYNCREALj),0); 





fillchar(SYNCIMAG, ,sizeof(SYNCIMAGn ),0); 
SYNCREALa[ 68] := -SYNCMAG; SYNCIMAG,{ 68]: = -SYNCMAG; 
SYNCREALa[ 69]: = -SYNCMAG; SYNCIMAG,[ 69]:= -SYNCMAG; 
SYNCREALa[ 70]: = -SYNCMAG; SYNCIMAGa[70]:= SYNCMAG; 
SYNCREAL«[ 71]:= -SYNCMAG; SYNCIMAG,[71]:= SYNCMAG; 
SYNCREALs[ 72]:= SYNCMAG; SYNCIMAGa,[ 72]:= -SYNCMAG; 
SYNCREAL,[ 73]:= SYNCMAG; SYNCIMAG,[73]:= SYNCMAG; 
SYNCREAL,[ 74]: = -SYNCMAG; SYNCIMAG,[{ 74]:= SYNCMAG; 
SYNCREALA[ 75]:= SYNCMAG; SYNCIMAGa[75]:= SYNCMAG; 
SYNCREALa[ 76]: = -SYNCMAG; SYNCIMAGa[76]:= SYNCMAG; 
SYNCREALA[ 77]: = -SYNCMAG; SYNCIMAG,[ 77]:= -SYNCMAG; 
SYNCREALa[ 78]: = SYNCMAG; SYNCIMAG~[78]:= -SYNCMAG; 
SYNCREALa[ 79]: = -SYNCMAG; SYNCIMAGa{ 79]:= -SYNCMAG; 
SYNCREALa{ 80]:= SYNCMAG; SYNCIMAGsf[ 80]:= SYNCMAG; 
SYNCREALA[ 81]: = SYNCMAG; SYNCIMAG~{ 81]:= -SYNCMAG; 
SYNCREALa[ 82]: = -SYNCMAG; SYNCIMAGa~{ 82]:= -SYNCMAG; 
SYNCREALA[ 83]: = SYNCMAG; SYNCIMAG,[{ 83]:= SYNCMAG; 
(*complex conjugate image) 
for J 3= 685to Gs do 
begin 


SYNCREALa[ J] ; 
-SYNCIMAGa[ J] ; 


SYNGREAMAZ sO |: 
SYNCIMAGa{ 256-J]: 
end; (*for J*) 


Comp lexFFT( 256, INVERSE, SYNCREAL,SYNCIMAG, ERROR); 


(*scale data/load time sequence*) 
for J := 0 to 255 do 
begin 
TEMP: =round(SYNCREALa[ J] + 126); 
if TEMP < 0 then 
TEMP: =0; 
SYNCDATA: =TEMP; 
BCST[ J+F IRST_ELEMENT] : =SYNCDATA; 
end; (*for J*) 


dispose(SYNCREAL); 

dispose(SYNCIMAG); 
end; (*SyncBaud*) 
( Moo oem OO OCe SO Oe Cee See Se SOS SHS Se SSS Se SSS Soe S ees Se SS SSS SSeS Se See Ses Seeeseeeoes we ) 
( ee ee ee * ) 


procedure SelectBaud; 
(*SelectBaud establishes kx, k1, and k2*) 


var 
ANSWER : integer; 


begin 
kx: =0; 
(*select baud size*) 
repeat 
if kx < 0 then writeln('TRY AGAIN'); 
writeln('What is the length of the bauds (kx)?'); 


Weitein( aoe. 256, 512,81024, 2048, 4096"); 
readin( ANSWER); 
case ANSWER of 
25 Ox 25105 
512: kx: =512; 
1024; kx:=1024; 
2048: kx: =2048; 
4096: kx: =4096; 
end; (*case ANSWER*) 
if kx = 0 then kx := -1; 
until kx > 0; 


(*set tone limits*) 
case kx of 
256: begin 
k1:=68; k2:=83; 
end; 
512: begin 
k1:=135; k2:=166; 
end; 
1024: begin 
k1:=269; k2: =332; 
end; 
2048: begin 
k1:=537; k2:=664; 
end; 
4096: begin 
role OV/ereeke2: — ls283 
end; 
end; (*case kx*) 
end; (*SelectBaud*) 
( Heme moeeeeocoo eee ee eeseocese ses eoe sss tesco sssoco esc eco feces eee sesso oesscecess bY g ) 


procedure TailorPacket; 
(*TailorPacket sets the maximum number of baud required to encode the 
message*) 


begin 
MESSAGESIZE: = filesize(BYTEFILE); 
writeln('Message is ',MESSAGESIZE,' bytes.'); 


(*kx/2 is the number of bit/baud for QPSK. kx/2-2 is the number 
for DQPSK. Each character is 8 bits*) 
CHARACTERS_PER_BAUD: =(kx/8 - 2)/8; 


(*61440/kx is maximum number of bauds possible*) 
MAXNUMCHAR: = trunc(61440. 0/kx * CHARACTERS_PER_BAUD); 
if MESSAGESIZE > MAXNUMCHAR then 
begin 
writeln('Message is to large. The last ', 
MESSAGESIZE - MAXNUMCHAR, 
characters will not be transmitted. '); 
MESSAGESIZE: =MAXNUMCHAR; 
end; 


MAXNUMBAUDS: =trunc(MESSAGESIZE / CHARACTERS_PER_BAUD); 
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(*ensure last few characters are included*) 
if frac(MESSAGESIZE / CHARACTERS_PER_BAUD) > 0.0 then 
MAXNUMBAUDS: =MAXNUMBAUDS + 1; 


repeat 

writeln; 

writeln('Enter number of ',kx,' bauds to process. ', 
MAXNUMBAUDS,' is the maximum. "); 

readln(NUMBAUDS); 


until NUMBAUDS in [1..MAXNUMBAUDS] ; 
end; (*TailorPacket*) 


procedure DiffEncode; 

(*DiffEncode differential encodes symbols on a tone-to-tone basis. 
BYTEFILE is read from one byte at a time. The byte is isolated into 
2-bit groups and stored in BITS. BITS is then used to DQPSK encode the 
frequency domain arrays XREAL and XIMAG. Bytes partially encoded are 
carried over into the next baud by global variable TEMPBYTE. *) 


var 
lf : integer; 
BITS : byte; 


begin 
fillchar(XREAL» ,sizeof(XREALn») ,0); 
fillchar(XIMAG~ ,sizeof(XIMAGa~) ,0); 


(*first tone of every baud set to pi/2*) 
XREAL«[ k1] : = MAGNITUDE; 
XIMAG~a[k1]: = MAGNITUDE; 


if SYMBOLCOUNT = 0 then 
Pedad( SYPEF ILE, VEMPRYTE 


(*break apart character”) 
for J:= (kl + 1) to k2 do 
begin 
SYMBOLCOUNT: =SYMBOLCOUNT + 1; 
if frac(SYMBOLCOUNT / 4) = 0.25 then 
BITS: = (TEMPBYTE and $CO) shr 6; 
if frac(SYMBOLCOUNT / 4) = 0.5 then 
BITS: = (TEMPBYTE and $30) shr 4; 
if frac(SYMBOLCOUNT / 4) = 0.75 then 
BITS: = (TEMPBYTE and $0C) shr 2; 
if frac(SYMBOLCOUNT / 4) = 0.0 then 
begin 
BITS: = TEMPBYTE and $03; 
if not EOF(CBYTEFILE) then 
read( BYTEFILE, TEMPBYTE) 
else 
TEMPBYTE: =$40; (*£i11 character*) 
end; 
if (BITS < 0) and (BITS > 3) then 
writeln('Bits not assigned properly'); 
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(*differential encode*) 
case BITS of 


2: begin XREALa[ J]:= XIMAGa[J-1]; 
XIMAGa[ J]: =-XREALaA[ J-1]; end; 
3: begin XREALa[ J] : =-XREAL«[ J-1]; 
XIMAGa[ J] ; =-XIMAGa[J-1]; end; 
end; (*case BITS*) 
end; (*for J*) 


O: begin XREAL,[J]:= XREALa[ J-1]; 
XIMAGa[J]:= XIMAGa[J-1]; end; 

1: begin XREALa[ J]: =-XIMAGa[ J-1]; 
XIMAGa[ J]:= XREALA[ J-1]; end; 


(*complex conjugate image*) 
for J:= kl to k2 do 
begin 
XREALs[ kx + J]:= XREALA[ J]; 
XIMAGa[ kx - J]: =-XIMAGa[ J]; 
writein(TESTFILE , BAUDCOUNT: 4,J:5,trunc(XREAL,[ J] ): 6, 
trunc(XIMAG,[ J] ): 6); 
end; 
end; (*DiffEncode*) 
( ie ee ee ee a eee ene ee ee ee ee on } 


procedure ScaleData; 

(*ScaleData converts each real value in XREAL down to a byte and stores 
the byte in the packet storage buffer, BCST. INDEX establishes the 
location in the buffer of each byte in the packet.) 


var 
IOSD oll, Ie ule : integer; 
DATA : byte; 
begin 
for J := 0 to kx-1 do 


begin 
TEMP round xREAls| J) + i2e); 
aise INSU! << (0) jelaveya 
TEMPS = m0. 
DATA := TEMP; 
(*256 is added to INDEX to start message bauds 
after the sync baud*) 
INDEX := J+(BAUDCOUNT-1)*kx+FIRST_ELEMENT+256; 
BCST[ INDEX] := DATA; 
end; (*for J*) 
end; (*ScaleData*) 


procedure Dmainit(var BCST : BCSTARRAY; 
BYTECOUNT : integer); 
(*Assembly language procedure used to initialize and unmask the DMA for 
data transfer. The source code must be converted to a BIN file.*) 


external 'DMAINITS. BIN'; 


procedure Dmastop; 
(*Masks DMA, stopping data transfer. *) 
external '‘DMASTOP. BIN’; 


(Wee nn nr rr rr rr re nr ener eee nen n ne *) 
( eee ee eee ee % ») 
begin 

clrscr; 


(*contains hex values to be encoded and transmitted”) 
assign( BYTEFILE, ‘MESSAGE. DAT'); 
reset(BYTEFILE); 


(*Output file of encoded symbols. Used for system testing*) 
assien@lesls lb, XMITDATADAN. ): 
rewrite(TESTFILE); 


INVERSE: =true; 
repeat 
writeln('Enter magnitude of tones.(must be > 0.0) '); 
readIn( MAGNITUDE); 
until MAGNITUDE > 0.0; 


writeln('Loading sync baud. '); 
SyncBaud; 


SelectBaud; 
TailorPacket; 


SYMBOLCOUNT: =0; 

TEMPBYTE: =S00; 

writeln('Number of bauds is ',numbauds); 
new(XREAL); 

new( XIMAG); 


for baudcount := 1 to numbauds do 

begin 
Dif fEncode; 
writeln('Performing IFFT ' ,BAUDCOUNT, ' a 

NUMBAUDS-BAUDCOUNT,' left’); 

Comp lexFFT( kx, INVERSE , XREAL, XIMAG , ERROR); 
ScaleData; 

end; (*for BAUDCOUNT*) 


dispose( XREAL); 
dispose(XIMAG); 
close(BYTEFILE); 
close(TESTFILE); 
BYTECOUNT := 256 + NUMBAUDS*kx - 1; 


repeat 
writeln('Press return to transmit’); read1n; 
Dmainit( BCST,BYTECOUNT); 
repeat 
writeln('Transmit some more? (*yes or no*) '); 
readln( ANSWER); 
until ANSWAR in ['n' jai, y 5 Yoel: 
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dmastop; 
until ANSWER in ['n','N']; 


(* reset(TESTFILE); 
while not EOF(TESTFILE) do 
begin 
while not EOLN(TESTFILE) do 
begin 
read( TESTFILE , NEXTCHAR); 
write(NEXTCHAR); 
end; 
readIn( TESTFILE); 
writeln; 
end; 
close(TESTFILE); *) 
end. 


APPENDIX G. DMAINIT 


codeseg segment 
public dmainit 
assume cs: codeseg 


; procedure DMAINIT ( BLKADDRESS : XMITPOINTER; 
‘ BYTECOUNT : INTEGER); 


;this procedure initializes dma channel 3 and sets the 

; parameters to output the array bcst by passing the address 
;of the array start on the stack. BYTECOUNT is the number 
;of bytes to transfer and is pushed on the stack by the 
;calling program 


dma equ 0 
dmapage equ 80h 
dmainit proc near 
push bp 
mov bp,sp ;use bp to address stack 
les di,dword ptr[ bpt+6];move add of best into es:di 
mov al,5bh ;dma chan 3 single mode, read, autoinit 
out dmat+ll,al 
out dmati2,al ;reset first/last ff 
mov ax,es ;calc high order 4 bits of buffer area 


mov cl,4 
rol ax,cl 


push ax 3save ax for dma start addr 

and al,Ofh 

out dmapaget2,al 3;store in ch 3 dma page reg 

pop ax 

and al,Of0h 

add ax,dd ;get page offset 

out dmat6,al ;output waveform buffer start addr 


mov al,ah 

out dmat6,al 

mov ax,[ bp+4] ;output dma byte count 
out dma+7,al 

mov al,ah 

out dmat+7,al 


mov a1,3 ;unmask ch 3 to start 

out dmat+10,al 

pop bp 

ret 6 ;pop 6 bytes off stack for addr of best 


dmainit endp 
codeseg ends 
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APPENDIX H. DMASTOP 
codeseg segment 
public dmastop 
assume cs: codeseg 
; procedure DMASTOP; 


3 
; this procedure stops dma channel 3 


dma equ 0 
dmapage equ 80h 
dmastop proc near 
push bp 
mov Alls ;mask ch 3 to stop 
out dma+10,al 
pop bp 
ret 


dmastop endp 
codeseg ends 


De 


APPENDIX I. RECEIVE 


program RECEIVE; 

(*This program acquires an analog signal, stores the raw data in a memory 
buffer, converts the data to TYPE real, performs an FFT using COMPFFT, 
differential decodes between adjacent tones and displays the result in 
the frequency domain with the phase quadrant represented by color. 

Data input can also be read in from the time domain sequence file 
HOUTFA1. DAT generated by program TRANSMIT) 


uses Graph, Crt, tp4d16; 


*#Sl<%* 


) 
(*$R-*) 


const 
kx = 256; (*set baudlength*) 
kxml = 255; 


type 
TNvector = array[{[0..kxml] of real; (*TYPE for real and imaginary 
data for FFT routing*) 
TNvectorPtr = aTNvector; (*Pointer for FFT data array 
which allows dynamic allocation 
of memory“) 


var 
INVERSE : boolean; 
XREAL, XIMAG : TNvectorPtr; 
ERROR : byte; 
NUMPTS, 
ki ie2 : integer; 
ANSWER : char; 


(“si EPRiS7B2, ING) 
(*S$T COMPRFT. ING) 


procedure AcquireData; 

(*AcquireData intializes Metrobyte DASH -16F data acquisition board. 
Using TTOOLS procedures D16_init and D16_ainm. Data transfer is 
controlled by the DMA controller and initialized by D16_ainm, and 
disabled by D16_dma_int_disable. TTOOLS procedures are external 
procedures included by ‘uses’ tp4d16.*) 


const max_buffer = 1000; 


var i: integer; 
rate: real; 
cnt_num, mode, cycle, trigger, 
base_adr, err_code, int_level, dma_level, 
board_num, chanlo, 
op_type, status, next_cnt, err_code_s : integer; 
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dmaPointer: pointer; 
datavector: ainteger; 
ad_data, chan_data: array[{0..max_buffer] of integer; 


begin 
clrscr; 


(*allocates memory on the heap and returns a pointer to the start 
of the buffer*) 
Get DMABuf fer(max_buffer,dmaPointer,err_code); 


(*This statement assigns a generic pointer to a variable of a 
specific pointer type, i.e. ainteger, so that the pointer can be 
passed to the dl6_ainm routine. *) 

datavector := dmaPointer; 

board_num := 0; int_level := 7; dma_level := 1; 

base_adr := $300; 


(*initialize the driver*) 

d16_init( board_num,base_adr,int_level,dma_level,err_code); 
chanlo := 0; 

cycle: =0; (*0 - one sweep of the DMA 1 - autoinitialize*) 
trigger: =0; (*0 - external 1 - internal) 

cnt_num: =kx; (*# of samples*) 

rate := 10000.0; (*used for internal trigger*) 

mode := 2; (*DMA mode*) 


writeln('Ready to acquire’); 

(*colects kx analog values using DMA and stores in a buffer*) 

dl6_ainm(board_num,chanlo,mode,cycle,trigger,cnt_num,rate, 
datavectora,err_code); 


(*status indicates the progress of acquisition. When all 
samples have been acquired status=0*) 
status := 11; 


(“wait until all data acquired*) 
repeat 
dl6é_dma_int_status(board_num,op_type,status,next_cnt, err_code_s); 
until status = 0; 
writeln('Data received'); 


if err_code <> 0 then 
d1l6_print_error(err_code) 
else 
begin 
writeln( ‘Processing data'); 


(*converts left justified data, returns the true binary value 
of the sampled data*) 

d1l6_convert_data(2047,cnt_num,datavectorn,ad_data[ 0], 

chan_data[ 0] ,0,err_code); 

new( XREAL); 

new( XIMAG); 

fillchar(XREAL,,sizeof (XREALj),0); 

fillchar(XIMAG~, ,sizeof (XIMAG~),0); 
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for I:= 0 to kxml do 
begin 
XREALA[ I] := ad_data[ I] /5; 
end; 
end; 


(*stop DMA operation if in autoinitialization*) 
di6_dma_int_disable(board_num,err_code); 


(*frees memory allocated with GetDMABuffer*) 
FreeDMABuffer(max_buffer,dmaPointer,err_code); 


end; (*Acquire*) 

(Me ee rn nn nee en nee n enna *) 
(aww mew ew rn nw ew er rn ee ne re ne een een ene ne een en nnn *) 
procedure GetData; 

(*GetData loads XREAL with the real time domain sequence generated by 
program TRANSMIT. Only one baud is loaded*) 


var 
TNS RR IS JOLIE. ‘tile of byte; 
DATA : byte; 
I, 
HEME : integer; 
begin 
assign(THEFILE, 'HOUTFA1. DAT' ); 
reset(THEFILE),; 
new( XREAL); 
new(XIMAG); 


fillchar( XREAL;» ,sizeof (XREALjA),0); 
fillchar(XIMAGs,sizeof (XIMAG~),0); 
for I := 0 to kxml do 


begin 
read(THEFILE, DATA); (*read file one byte at a time*) 
TEMP := DATA; (*puts byte into integer variable”) 
XREALa[I] := (TEMP - 126); (*loads XREAL array*) 
end; (*for I*) 
close(THEFILE); 
end; (*GetData*) 
(Wea nn a en ne nen nn nn nee eee eee eee ee eee ee eee sees *) 
(ee * 


procedure DiffDecode; 

(*DiffDecode differentially decodes complex frequency domain arrays XREAL 
and XIMAG. Four decoded symbols are recombined into a byte and 
transferred to file BYTESOUT. DAT. *) 


var 
I, 
SYMBOLCOUNT : integer; 
TEMPREAL,TEMPIMAG : real; 
ICIS) 5 Mlepaledes ee. : byte; 
TEMPCHAR 9 lausicx 
OUTFILE elles 
begin 
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assign(OUTFILE, 'BYTESOUT. DAT'); 
rewrite(OUTFILE); 
SYMBOLCOUNT: =0; 
TENE TE? 0: 
for 1 := ki to (k2 - 1) do 
begin 
(*Complex multiply two adjacent tones, I and the complex conjugate of 
I+l. This will give the phase difference between the two tones. 
The answer is in rectangular notation*) 
TEMPREAL: =XREALa[ I] * XREALa[I+1] + XIMAGa[I] * XIMAG~[I+1]; 
TEMPIMAG: =XREALa[ I] * XIMAGa[I+1] - XREAL«[I+1]) * XIMAG~[I]; 


(*Complex multiply (TEMPREAL + j TEMPIMAGE) and (1+j). This rotates 
the differential vector 45 degrees. XREAL and XIMAG are used to 
store the results. This eliminate the original data*) 

XREALa«[ I] : =(TEMPREAL - TEMPIMAG)/80; 
XIMAGa[ I] :=(TEMPREAL + TEMPIMAG) /80; 


if (XREAL«~[I] >= 0) and (XIMAG,[I] > 0) then 
BITS: =$00; 

if (XREAL«[I] < 0) and (XIMAGa[I] > 0) then 
BITS: =$01; 

if (XREAL-[I] < 0) and (XIMAGa[I] <= 0) then 
BITS: =$03; 

if (XREAL«[I] >= 0) and (XIMAG-[I] <= 0) then 
BITS: =$02; 

SYMBOLCOUNT := SYMBOLCOUNT + 1; 


(*fill TEMPBYTE with four symbols*) 
if frac(SYMBOLCOUNT / 4) = 0.25 then 
TEMPBYTE := (BITS shi 6); 
if frac(SYMBOLCOUNT / 4) = 0.5 then 
TUEMIOEN INE, SE (CRBS Gleull )) pe Washers citia,s 
if frac(SYMBOLCOUNT / 4) = 0.75 then 
TUESUEBNTINE, SS (CRUCUS einil 2) Gre Welleshavee 
if (frac(SYMBOLCOUNT / 4) = 0.0) then 
begin 
TEMES, 2S LIN ere INS INES 
TEMPCHAR := chr( TEMPBYTE); 
write(OUTFILE,TEMPCHAR); 
TEMPBYTE: =0; 
end; (*if frac*) 
end; (“for I*) 


TEMPCHAR := chr( TEMPBYTE); 
write(OUTFILE,TEMPCHAR); (*puts two 0 at end of each baud*) 
close(OUTFILE); 


XREAL«[ k2] : =0; 

XIMAGaA[ k2] : =0; 
end; (*DiffDecode*) 
( New eww mm ewe een eee ee ww eee eee eee eee ee eee ee eee eee eer *) 
procedure GraphData(XR, XI: TNVectorPtr); 


(*GraphData graphs complex arrays XR and XI. Two scales are used: Full 
scale, 0 to kx/2 tones and Zoom scale, kl to k2*) 


63 


var 
grDriver, grMode,ErrCode, 
Tae (*loop indices*) 
kel kz, 
MXDIV2, 
COLOR, 
Ake, 
Woy ME : integer; 


begin 
ki :="round(kx * 67.0 / 256: 0): 
k2 := round(kx * 83.0 / 256.0); 
grDriver := EGA; 
grMode := EGAHi; 
initgraph( grDriver,grMode,'c: TP4 GRAPHICS'); 
ErrCode := GraphResult; 
if ErrCode <> grOK then 
begin 
writeln('Graphics error: ', GraphErrorMsg(ErrCode)); 
writeln('Program aborted...'); 
Halt(1); 
end; (*if ErrCode*) 
SetBKColor( black); 
SetColor(white); 


(*axis contruction*) 
line( 60,300,600, 300); (*x-axis for full spectrum’) 
line(68,140,600,140); (*x-axis for zoom spectrum”) 
line(332,300,332,308); (*x-axis scale marks*) 
line(392,300,392,308); (*x-axis scale marks*) 
line(60,300,60,308); 
line(572,300,572, 308); 
line(96,140,332,300); 
line(576, 140,392,300); 


MXDIV2 := round(kx / 2); 


for I := 0 to MXDIV2 do 
begin 
if (XRa[I] >= 0) and (XIa[I] > 0) then 
SetColor( lightmagenta); 
if (XRa[I] < 0) and (XI,[I] > 0) then 
SetColor(lightgreen); 
if (XRa[I] < 0) and (XI,[I] <= 0) then 
SetColor(lightred); 
if (XR»[I] >= 0) and (XIa[I] <= 0) then 
SetColor( yellow); 


(*zoom spectrum* ) 


if (I > kl) and (I <= k2) then 
begin 
XZ := round(2*((I * 4096.0 / kx) - 1040)); 
YZ := 140 - round(sqrt(sqr(XRa[I]) + sqr(XIa[I] ))); 
line( XZ,140,XZ2,YZ); 


end; 
(*full spectrum*) 


X := round(2*(I1/MXDIV2 * 256 + 30)); 
Y := 300 - round(sqrt(sqr(XRa[I]) + sqr(XIa[I]))); 
line(X,300,X,Y); 

end; 


(*graph labels*) 
setcolor(lightgray), 
outtextXY(60,320,'0'); 

Outtextny¢osl 320, ki’ ys 
outtextXY(391,320,'k2'); 

outtextXY(565 ,320,'kx/2'); 

outtextXY(150,330,' frequency--->'); 

outtextXY( 20,250, amplitude'); 

readln; 

outtextXY(300,330,'Press return to continue. '); 


readln; 
CloseGraph; 
end; (*GraphData*) 
(Bann nnn nn nn nn nn nn ne nn nnn een eee e ncn e en neeeee *) 
(* a eae *) 
begin (*main body*) 
INVERSE := false; (*sets forward FFI*) 
ERROR := 0; 
NUMPTS := kx; 
kil S= seemunel(lee “ C750 // 2550 se ile 
k2 := round(kx * 83.0 / 256.0); 


repeat 
writeln('Enter 1 to sample data’); 
writeln(' 2 to read data from an ASCII file'); 
readin( ANSWER); 

DicreANSWER am { 1°, 2°]; 


if ANSWER = '1' then 
(*AcquireData samples input analog signal*) 
AcquireData 
else 
GetData; 

(*GetData reads a file of discrete sample produced in the transmitter. 
GetData is used for testing Graphics and Decoding procedures without 
the need for external hardware and Data aquisition board 
initialization*) 


writeln('Computing FFT'); 
ComplexFFT (NUMPTS, INVERSE, XREAL, XIMAG, ERROR); 
DiffDecode; 
(* writeln('Error = ',ERROR,' hit the enter key');readln;*) 
GraphData(XREAL, XIMAG); 
end. 
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APPENDIX J. RECMES 


program RECEIVER; 

(“Acquires the signal. Stores it in a memory buffer. Differential decodes 
between tones. Maximum number of bauds are received. The number of 

bauds processed is a user input*) 


uses Graph, Crt, tp4d16; 


(*SI-*) 
(*SR-*) 


const MAX_BUFFER = 65500; 


type 
TNvector = array-0..4096- of real; (“TYPE for real and imaginary data 
for FFT routing”) 
TNvectorPtr = aTNvector; (*Pointer for FFT data array which 
allows dynamic allocation of memory”) 


var 
INVERSE : boolean; 
XREAL, XIMAG :TNvectorPtr; 
ERROR, TEMPEYIDE byte; 
J, 
k1,k2,kx,ANSWER, 
ERR_ COLE, “BAUDCOUNT ; 


SYMBOLCOUNT , 

NUMBAUDS , 

MAXNUMBAUDS : integer; 
MAGNITUDE, 

PHASE : real; 
DATAVECTOR > ainteger; 
DMAPOINTER : pointer; 
TESTE ILE, 

OUTFILE rE; 


C*Sl Pilar be ence.) 
(*$1I COMPFFT. INC*) 


(* we a a ww a a a a a a ee ee ee we ew ew ee ee ee ee ee wee ee ee ee ea ee ee He eee ee eee ee we) 
(* we ww www www wwe se ee ee www www wwe eee et wee ww em eee ewww eee eee eee eee ee eee eee ee *) 
procedure PacketSetUp; 
begin 
repeat 
Clissiea; 


if kx < 0 then writeln('TRY AGAIN'); 
writeln(' Enter baud size '); 
readln( ANSWER); 
case ANSWER of 
25 6 exes 250; 
Siz: hepa S12; 
1024: kx: =1024; 
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2048: kx: =2048; 
4096: kx: =4096; 
end; (**case*) 
if kx = 0 then kx := -1,; 
Misia, tee S Oe 
MAXNUMBAUDS := trunc( (MAX_BUFFER/2)/kx); 


repeat 
writeln; 
writeta( Enter number of kx, bauds to process. , 
MAXNUMBAUDS,' is the maximum. '); 
readln( NUMBAUDS); 
until NUMBAUDS in [1..MAXNUMBAUDS] ; 


kl := round(kx * 67.0 / 256.0 + 1); 


k2 := round(kx * 83.0 / 256.0) 
end; (*PacketSetUp*) 
(Hala cio omens ec ew een ee ew wa ean aa aero ewe nea n ea sen sseeaseeserns=-e=--= *) 
( Ce ee ) 


procedure AcquireData; 

(*AcquireData initializes Metrobyte DASH-16F data acquisition board, 
using TTOOLS procedure Dl6_int and Dl6_ainm. Data transfer is 
controlled by the DMA controller and initialized by D1l6_ainm and 
disabled by D16_dma_int_disable. TTOOLS procedures are external 
procedures included by ‘uses’ tp4d16.*) 


var 

RATE: real; 

I,CNT_NUM, MODE, CYCLE, TRIGGER, 

BASE_ADR, INT_LEVEL, DMA_LEVEL, 

BOARD_NUM, CHANLO, 

OEE YEE RS TARU SME NE Mia NiaEer Ra OD Ems. integer; 
begin 


BOARD_NUM := 0; INT_LEVEL := 7; DMA_LEVEL := 1; 
BASE_ADR := $300; 


D16_init(BOARD_NUM, BASE_ADR, INT_LEVEL,DMA_LEVEL,ERR_CODE); 


CHANLO := 0; 

CYCLE: =0; (*O-one sweep of the DMA 1l-autoinitialize*) 
TRIGGER: =0; (*0 - external 1 - internal*) 

CNT_NUM: =trunc(MAX_BUFFER / 2); (*# of samples*) 

RATE := 10000. 0;(*used for internal trigger*) 

MODEM = 25 (*DMA mode*) 

writeln('Ready to acquire’); 


D16_ainm( BOARD_NUM , CHANLO , MODE , CYCLE , TRIGGER, CNT_NUM, 
RATE, DATAVECTOR»~ , ERR_CODE); 


STATUS := 11; 
(*status indicates the progress of acquisition. When all 


samples have been acquired status=0*) 


repeat 
D16_dma_int_status(BOARD_NUM,OP_TYPE ,STATUS ,NEXT_CNT, ERR_CODE_S); 
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until STATUS = 0; 
writeln('Data received’); 


if ERR_CODE <> 0 then 
D16_print_error(ERR_CODE); 
D16_dma_int_disable( BOARD_NUM,ERR_CODE); 


end; (*Acquire*) 

Ree eee ee ee nw ee ee ee eee ne ee ee ee ee ew ee ee ee ee ee ee ee ee ee ee ee * 
procedure ConvertData; 

(*ConvertData seperates channel and acquired data. CHAN_DATA is not used. 
Acquired data is stored in XREAL.*) 


var 
AD_DATA: array[0..4096] of integer; 
I ,;CHAN_DATA ,ERR_CODE, 
SEGMENTPART ,OFFSETPART: integer; 
NEWDATAVECTOR :ainteger; 
TEMPPOINTER : pointer; 
begin 


writeln('Processing baud ' ,BAUDCOUNT); 

fillchar( XREAL;»,sizeof (XREAL~s),0); 

fillchar(XIMAG~s ,sizeof (XIMAG~),0); 

SEGMENTPART: =seg(DATAVECTOR,~ ); 

OFFSETPART: =ofs(DATAVECTOR,) + 2 * kx * (BAUDCOUNT - 1); 

TEMPPOINTER: =ptr( SEGMENTPART , OFFSETPART); 

NEWDATAVECTOR := TEMPPOINTER; 

d16_convert_data( 2047 ,kx ,NEWDATAVECTOR,«~ ,AD_DATA[ 0] , 
CHAN_DATA,0,ERR_CODE); 


for 1:= 0) to (kx = 1) do 
begin 
XREALa[ I] := AD_DATA[ I] /5; 
end; 
end; (*ConvertData*) 


procedure DiffDecode; 

(*DiffDecode differentially decodes complex frequency domain arrays XREAL 
and XIMAG. Four decoded symbols are recombined into a byte and 
transferred to file BYTESOUT. DAT. *) 


var 
I : integer; 
TEMPREAL,TEMPIMAG : real; 
BITS : byte; 
TEMPCHAR : char; 
begin 


for I := kl to (k2 - 1) do 
begin 
(*Complex multiply two adjacent tones, I and the complex conjugate of I+1. 
This will give the phase difference between the two tones. The answer 
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is in rectangular notation*) 
TEMPREAL: =XREALa[ 1] * XREALA[ I+1] + XIMAGa[I] * XIMAGa[I+1]; 
TEMPIMAG: =XREALa[I] * XIMAGa[I+1] - XREALs[I+1] * XIMAG~[ I]; 


(*Complex multiply (TEMPREAL + j TEMPIMAGE) and (1+j). This rotates the 
differential vector pi/4 radians. XREAL and XIMAG are used to store 
the results. This eliminate the original data*) 

XREAL«[ I]: =(TEMPREAL - TEMPIMAG)/80; 
XIMAGa[ I]: =(TEMPREAL + TEMPIMAG)/80; 


(*decode* ) 

if (XREAL»[I] >= 0) and (XIMAGa[I] > 0) then 
BITS: =$00; 

if (XREALs[ I] < 0) and (XIMAGa[I] > 0) then 
BITS: =$01; 

if (XREALs[I] < 0) and (XIMAGa[I] <= 0) then 
BITS: =$03; 

if (XREALs[ I] >= 0) and (XIMAG,[I] <= 0) then 
BITS: =$02; 


SYMBOLCOUNT := SYMBOLCOUNT + 1, 


(*fi11 TEMPBYTE with four symbols*) 
if frac(SYMBOLCOUNT / 4) = 0.25 then 
EMER Ee =" (CBistsms le 6): 
if frac(SYMBOLCOUNT / 4) = 0.5 then 
TEMPBYTE := (BITS shl 4) or TEMPBYTE; 
if frac(SYMBOLCOUNT / 4) = 0.75 then 
TEMPBYTE := (BITS shl 2) or TEMPBYTE; 
if (frac(SYMBOLCOUNT / 4) = 0.0) then 
begin 
TEMPBYTE := BITS or TEMPBYTE; 
TEMPCHAR := chr( TEMPBYTE); 
write(OUTFILE , TEMPCHAR); 
TEMPBYTE: =0; 
end; (*if frac*) 
end; (*for I*) 


XREAL»[ k2]: 
XIMAGa[ k2] : 


=i, (08 
=1.0 


. 
a) 


end; (*DiffDecode*) 
ee er rr roe * ) 
procedure Showmessage; 

(*Showmessage reads in decoded message*) 


var 
NEXTCHAR: char; 


begin 

writeln; 
writeln('The message transmitted is..'); 
assign(OUTFILE, ‘MESSAGE. DAT" ); 
reset(OQUTFILE); 
while not EOF(QUTFILE) do 

begin 

while not EOLN(OUTFILE) do 
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begin 
read( OUTFILE , NEXTCHAR); 
write(NEXTCHAR); 
end; (*while not EOLN*) 
readlin( OUTFILE); 
writeln; 
end; (*while not EOF*) 
close(OUTFILE); 
end; (*Showmessage*) 
(* we a ww os a a a a ww ww ww ww we ww ew ew we ee we ee ee ee ee ee eee ee eee eee eee *) 


begin (*main body*) 
GetDMABuf fer(MAX_BUFFER, DMAPOINTER ,ERR_CODE); 


DATAVECTOR := DMAPOINTER; (*This statement assigns a generic pointer 
to a variable of a specific pointer type, 
i.e. ainteger, so that the pointer can be 
passed to the d1l6_ainm routine. *) 


new( XREAL); 
new(XIMAG); 


INVERSE := false; (*sets forward FFT*) 


ERROR := 0; 
kx: =0; 


PacketSetUp; 


(*received message output file”) 
assign(OUTFILE, 'MESSAGE. DAT'); 
rewrite(OUTFILE); 


(*RECSTAT file shows the computed real and imaginary values of the 
received signal and its decoded representation*) 

assign(TESTFILE, 'RECSTAT. DAT"); 

rewrite(TESTFILE); 
(*  writeln(TESTFILE,'Baudlength is',kx:5);*) 


SYMBOLCOUNT: =0; 
TEMPBYTE: =0; 


AcquireData; (*AcquireData samples input analog signal*) 


for BAUDCOUNT := 1 to NUMBAUDS do 
begin 
ConvertData; 
Comp lexFFT (kx, INVERSE , KREAL, XIMAG ,ERROR); 


(* writeln(TESTFILE, 'Baud': 11, 'Tone :5,'Real’: 10, 
"Imaginary': 10,'Mag': 10, Phase': 8); *) 
EO elo eet om ecm e 
begin 
MAGNITUDE: =20* 1In(sqrt(sqr( XREALa[ J] )+ 
sqr(XIMAGa[J]))) / 1n(€10); 


if (XREALa[J] >= 0) then 
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PHASE: =arctan(XIMAGa[ J] / XREAL«[ J] ); 

if (XREALs[J] < 0) and (XIMAGs[J] > 0) then 
PHASE: =Pi + arctan(XIMAGa[ J] / XREALa[J] ); 

if (XREAL»[J] < 0) and (XIMAG,[J] <= 0) then 
PHASE: =arctan(XIMAGA[J] / XREALs[J]) - Pi; 


writeln(TESTFILE,'R' ,BAUDCOUNT:4,J:5, 
XREALa[ J]: 10: 4,XIMAGa[ J]: 10:4, 
MAGNITUDE: 10: 4, PHASE: 8: 4); 
end; 


DiffDecode; 


for J:=kl1 to k2 do 
begin 
MAGNITUDE: =sqrt(sqr(XREALa[ J] )+sqr(XIMAGa[ J] ))5 
write in(TESIFILE, 0 ,BAUDCOUNT: 4 ,J: 5, 
XREAL«[ J] /MAGNITUDE: 10:4, 


XIMAGa[ J] /MAGNITUDE: 10: 4); 
end; 
end; 


writeln(OUTFILE,chr( TEMPBYTE)); 


close(OUTFILE); 

close( TESTFILE); 

dispose( XREAL); 

dispose( XIMAG); 

FreeDMABuf fer(MAX_BUFFER , DMAPOINTER,ERR_CODE); 
Showmes sage; 


(* writeln( ‘Error = ',ERROR,' hit the enter key’); read1n;*) 
end. 


71 


APPENDIX K. RESPONSE 


program RESPONSE; 

(*RESPONSE outputs the channels phase response, by subtracting 
the transmitted phase from the received phase for each tone. 
Prior to running, set NUMTONES and K1MINUS1 for the proper 
baud parameters then compile and run the program*) 


var 
MAG, PHASE, 
XREAL,XIMAG : array[1.. 256] of real; 
K1MINUS1, 
NUMTONES ,K, 
I,BAUD,Kx,Kr : integer; 
XMITPHASE, 
XR,XI g e@Glg 
CODE : char; 
XMITDATA, 
RECDATA, 
CHANNELDATA : text; 
begin 
assign(XMITDATA, 'xmitdat1. 024'); 
reset(XMITDATA); 


assign(RECDATA, 'recstat.dat'); 
reset(RECDATA); 


assign(CHANNELDATA, 'resp1024. dat'); 
rewrite(CHANNELDATA); 


(*NUMTONES and K1MINUS1 should be set for the same baud parameters 
as xmitdat*. *** and recstat. dat*) 
NUMTONES: =64; 
K1MINUS1: =268; 


(*get transmitted phase and received phase*) 
while not EOF(XMITDATA)do 
begin 
for I:=1 to NUMTONES do 
begin 
readln( XMITDATA, BAUD ,Kx ,XREAL[ 1] ,XIMAG[ i} ); 
readln(RECDATA, CODE, BAUD, Kr,XR,XI,MAG[ I] , PHASE[ I] ); 
end; 


(*ensure tones match*) 
if Kx <> Kr then 
begin 
writeln('Kx <> Kr'); 
readin; 
end; 


(*convert transmitted symbol into phase*) 
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for I:=1 to NUMTONES do 
begin 

if (XREAL[I]>= 0.0) and (XIMAG[1I]>= 0.0) then 
XMITPHASE: = pi/4; 

if (XREAL[I]< 0.0) and (XIMAG[ i] >= 0.0) then 
XMITPHASE: = 3*pi/4; 

if (XREAL[ I] >= 0.0) and (XIMAG[I]< 0.0) then 
XMITPHASE: = -pi/4; 

if (XREAL[I])< 0.0) and (XIMAG[I]< 0.0) then 
XMITPHASE: = -3*pi/4; 


(*subtract transmitted phase from received phase*) 
PHASE[ I]: =PHASE[I] - XMITPHASE; 


(*phase differences may jump be 2*pi radian. This helps 
but does not work in every case. Output file may have 
to be edited to 2*pi jumps*) 

i if PHASE[I] > 2*pi then 
PHASE[ i]: =PHASE[I] - 2*pi; 
TEePHASE] i] < -2 then 
PHASE[ i]: =PHASE[ i] + 2*pi; 


(*actual tone numbers*) 
K:=I + K1MINUS1; 


writeln(CHANNELDATA, BAUD: 4,K: 4,MAG[ I]: 10: 4,PHASE[ I]: 10:4); 

readln(RECDATA); 

end; 
end; 

close( XMITDATA); 

close(RECDATA); 

close( CHANNELDATA); 


end. 
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APPENDIX L. STATISTICS 


program Statistics; 

(*to run this program the transmitter program Xmitmes must be run using 
MESSAGE. DAT. The file produced, XMITDAT.DAT, must be transferred to the 
receiver computer. Then transfer the data between computers using 
Xmitmes and Receive. Receive produces the file RECSTAT. DAT. This 
program uses XMITDAT. DAT and RECSTAT. DAT to compute the mean and 
variance for each quadrant, and the SNR for each quadrant, baud and 
the overall SNR for the packet.*) 


uses CRT; 


type 
STATARRAY = array[1..125] of real; 


var 
Kx ,K, NUMBAUDS , 
1,J,BAUD, 
TONEX , TONER, 
Q1,Q2,Q3,Q4, 
BITERROR1, BITERROR2, BITERROR3 , BITERROR4 : integer; 


TEMPR,TEMPI , 

SNRINDB, SNRIN, 

SNROUTDB, SNROUT, 

MEANR , MEANI, VARR, VARI: real; 


CODE : char; 
XD,RD,SD > text; 


KR, XIRR.RI 2array| ie. 16,127. 128] ot real 
BIM,BIV,BTQ : array{1..16] of real; 
Shigokee oh ogeka, 

SI1,S12,S13,S14 : STATARRAY; 


STATMATRIC : array[1..16,1..4,1..5] of real; 


Cee Oooo Sooo oe oonoo soe cen oop soo Se one so Sess sos soso eos e Hoss csccsosssoescs *) 


procedure Stat(NPTS: integer; 
X: STATARRAY; 
var XMEAN, XVAR : real); 


(“input NPTS:number of points to be used to compute the mean 
X: the array of data 


output XMEAN: the mean of the data in array X 
XVAR: the variance of the data in array X*) 


var 


SUM: real; 
N : integer; 
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begin 
(*Compute mean*) 
SUM: = 0.0; 


for N:= 1 to NPTS do 
SUM: = SUM + X[N}; 


if NPTS > O then 
XMEAN: = SUM / NPTS 
else 
SMEAN: =1. 0e-10; 


(*Compute variance*) 
SUM: = 0. 0; 


for N:=1 to NPTS do 
SUM: = SUM + sqr(X[N] - XMEAN); 


alte INDE IS) S> IL felneua 

XVAR: = SUM / (NPTS - 1) 
else 

XVAR: =1. 0e-10; 


end; (*Stat*) 


(? ae a ae ere ee a a ee en ee ee ee eal a wy 
begin (“main body*) 
repeat 
write('Enter baud size. '); 
readln( Kx); 


until (Kx mod 256 = 0) and (Kx >= 256) and (Kx <= 4096); 


write('Enter input SNR in DB '); 
readin(SNRINDB); 
SNRIN := exp(SNRINDB / 10 * 1n(10)); 


K := trunc(Kx / 16); 


assign(XD, 'XMITDAT. DAT'); 
reset(XD); 
assign(RD, 'RECSTAT. DAT'); 
reset(RD); 
case Kx of 


256: begin assign(SD,'STAT256. DAT'); NUMBAUDS:= 16; end; 
512: begin assign(SD,'STAT512. DAT'); NUMBAUDS:= 8; end; 
1024: begin assign(SD,'STAT1024. DAT' );NUMBAUDS:= 4; end; 
2048: begin assign(SD, 'STAT2048. DAT' ); NUMBAUDS: = Zens 
4096: begin assign(SD,'STAT4096. DAT');NUMBAUDS:= 1; end; 


end; (*case Kx*) 
rewrite(SD); 


(*read in transmitted and received data*) 
for I:= 1 to NUMBAUDS do 
begin 
for J:= 1 to K do 
readin(RD); 
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for J:= 1 to K do 
begin 
readin(XD, BAUD, TONEX, XR[I,J}, XI[I,J]); 
readin(RD, CODE, BAUD, TONER, RR[I,J], RI[I,J]); 
if TONEX <> TONER then 
begin 
clrscr; 
writeln(’XMITDAT. DAT and RECSTAT. DAT are in error. '); 
Halt; 
end; 
end; (*for J*) 
end; (*for I*) 
close(XD); 
close(RD); 


(*decode transmitted data*) 
for I:= 1 to NUMBAUDS do 
for J:= 1 to (K - 1) do 
begin . 
TEMPR: = XR[1,J] * MRE 1,J+1]) + Xu) Ae 
TEMPI:= XR[1,J] * XIL1,J+1] = ARPD +1)" = xia: 
XE 120] >= TEMPRS— TEMP a: 
XI[ 1,3) := TEMPR + TEME Td; 
end; (*for J*) 


(*store statistics for each quadrant and baud compute statistics*) 
BITERRORI: =0; BITERROR2: =0; BITERROR3: =0; BITERROR4: =0; 
for I:= 1 to NUMBAUDS do 
begin 
Qi: =0; Q2:=0; Q3:=0; Q4:=0; 
for J:= 1 to (K - 1) do 


begin 
Gf (XR[ 1,3) 2= 0.0) and (Xie 20..0)) then 
begin 
Ol: =Ol aie 
SRIPQUs = Relea) 
STIpei: = ipl, J]; 
if RE[ I,J) < 0.0 then 
BITERROR1: = BITERROR1 + 1; 
if RI[l J] <= 0.0 then 
BITERROR1: = BITERROR1 + 1; 
end; 
if (XR[ I,J] = 000) and CXI[i34]> 0.0) then 
begin 
Q2:=@2, + 15 
SR2(02]- = RR 1, J]: 
S12[@2)|: = RI 1,3); 
if RR[I,J] >= 0.0 then 
BITERROR2? =" BITERROR2 + 1; 
if RI{I,J] <= 0.0 then 
BITERROR2: = BITERROR2 + 1; 
end; 
if (XR[I,J] < 0.0) and (XI[1,J] <= 0.0) then 
begin 


Q3:=Q3 + 1; 
SR3[ Q3]:= RR{[I,J]; 
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SI3(Qsie= Ki 1.315 

if RR[I,J] >= 0.0 then 
BITERROR3: = BITERROR3 + 1; 

if Ki[ lag) > 0.0 then 
BITERROKS: = BITERRORGe+ 1; 


end; 
if (XR[I,J] >= 0.0) and (XI[I,J] <= 0.0) then 
begin 
Q4:=Q4 + 1; 
SR4[ Q4]:= RR[ I,J]; 
SI4[Q4]:= RI[I,J]; 
if RR[I,J] < 0.0 then 
BITERROR4: = BITERROR4 + 1; 
if RI[1I,J] > 0.0 then 
BITERROR4: = BITERROR4 + 1; 
end; 


> 
end; (*for J*) 


(*STATMATRIC is a 3 dimensional array. The first index is the baud 
number. The second is the quadrant. The third is the statistics 
for each quadrant. 
STATMATRIC[ I,J, 1] 
SRATMATRIG! LJ, 2] 
SHALTHATRIC| 153,53] 
STATMATRIC[ I,J, 4] 
STATMATRIC[ I,J, 5] 


the mean of the real parts 

the variance of the real parts 

the mean of the imaginary parts 

the variance of the imaginary parts 

the number of points transmitted in a given 
quadrant*) 


BGM 1) =020; 
BTV[ I]: =0. 0; 
BTQ[ I] : =0. 0; 
if Q1 > 1 then 
begin 
Stat(Q1,5R1,MEANR, VARR); 
Stat(Q1,SI1,MEANI, VARI); 
STATMATRIC[ I,1,1] : =MEANR; 
STATMATRIC[ 1,1,2]:=VARR; 
STATMATRIC[ I,1,3]: =MEANT; 
STATMATRIC[ 1I,1,4]:=VARI; 
BIM[ I]: =BIM[ I] + Q1*( abs(MEANR)+abs(MEANT)); 
BTV[ I]: =BTV[ I] + €Q1-1)*(VARR+VARI); 
BIOPRs=Bi0( 1] +22701; 
end, 
STATMATRIC[1,1,5] :=Q1; 


if Q2 > 1 then 
begin 

Stat(Q2,SR2,MEANR,VARR); 
Stat(Q2,S5I2,MEANI,VARI); 
STATMATRIC[1,2,1] : =MEANR; 
STATMATRIC[ 1,2, 2] : =VARR; 
STATMATRIC[ 1,2,3] : =MEANT; 
STATMATRIC[ 1,2,4] : =VARI; 
BTM[ I]: =BTM[ 1] + Q2*( abs(MEANR)+abs(MEANI)); 
BTIV[ I]: =BTV[ I] + (Q2-1)*(VARR+VARI); 
BrQ[a]: =BTQ[ I] + 2*@2; 
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end; 
STATMATRIC( 1, 2,5]: =@23 


if Q3 > 1 then 
begin 


Stat(Q3,SR3,MEANR, 
Stat(Q3,S13,MEANI, 


STATMATRIC[1,3,1]: 
STATMATRIC[ I,3,2]: 
STATMATRIC[ I,3,3]: 
STATMATRIC[ 1,3,4]: 


VARR); 
VARI); 
=MEANR; 
=VARR; 
=MEANT; 
=VARI; 


BTM[ I]: =BTM[1] + Q3*(abs(MEANR)+abs(MEANI)); 
BIV[ 1]:=BTV[I] + (Q3-1)*(VARR+VARI); 
BTQ[ 1]:=BTQ[I] + 2*Q3; 


end; 
STATMATRIC[ 1,3,5]: =Q3; 


if Q4 > 1 then 
begin 
Stat(Q4,SR4,MEANR,VARR); 
Stat(Q4,SI14,MEANI,VARI); 
STATMATRIC[ 1,4,1] : =MEANR; 
STATMATRIC[ 1,4, 2]: =VARR; 


STATMATRIC[ 1,4,3]:=MEANTI; 
STATMATRIC[ 1,4,4]:=VARI; 
BIM[ 1]: =BTM[ 1] + Q4*(abs(MEANR)+abs(MEANI)); 
BTV[ 1]: =BTV[I] + (Q4-1)*(VARR+VARI); 
BIQ[ 1]: =BIQ[ 1] + 2°04; 
end; 


STATMATRIC[ 1,4,5] : =Q4; 
end; (*for I*) 


(*output statistics to a file*) 


writeln(SD,'*** SNRIN= ',SNRIN: 10:7,' or ',SNRINDB:10:7,' DB 


writeln(SD); 
writeln(SD,' Baud length is ',Kx); 
SR2[ 1]: =0. 0; 
SR2[ 2]: =0. 0; 
SR2[ 3] : =0. 0; 
for I:= 1 to NUMBAUDS do 
begin 


writeln(SD); writeln(SD); writeln(SD);writeln(SD); 


sedese ) 


writeln(SD,'DFT statistics for baud #: ',1I:4);writeln(SD); 


LOGO — lomo 
begin 
(*real output data*) 
writeln(SD);writeln(SD); 


writeln(SD,'Given the phase is in quadrant : 


writeln(SD); 
if STATMATRIC[1,J,5] > 1 then 
begin 
writeln(SD,'Number of points = ':35, 
STATMATRIC[ 1,J,5]: 10:5); 
writeln(SD,'Mean of real parts = ':35, 


STATMATRIC[ I,J,1]: 10:5); 
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writeln(SD,'Variance of real parts = ':35, 
STATMATRIC[ 1I,J,2]: 10:5), 

SNROUT: =sqr(STATMATRIC[I,J,1]) / STATMATRIC[ I,J,2]; 

SNROUTDB:= 10.0 * In(SNROUT) / 1n( 10.0); 


writeln(SD,'SNROUT of real pact = 535, 
SNROUT: 10:5,' =", SNROUTDB: 10:5, “2nB ye 
write1n(SD, 'SNROUT/SNRIN Cheeccal part — 35, 


(SNROUT / SNRIN): 10:5); 


(*imaginary output data*) 

writeln(SD); 

writeln(SD,'Mean of imag parts = ':35, 
,STATMATRIC[T, 5]: 10: 535 

writeln(SD,'Variance of imag parts = ':35, 
STATMATRIC[I,J, we OS) 

SNROUT: = sqx(STATMATRIC[ I,J noe STATA Tait ,),4) ; 

SNROUTDB:= 10 * 1n(SNROUT) / In(10, 0); 


writeln(SD,'SNROUT of imag part = :35, 
SNROUT: 10:5,' = ', SNROUTDB: 10: Ba) Wee 
writeln(SD, ' SNROUT/SNRIN ef imae part = :35, 


(SNROUT / SNRIN): 1025); 
end (*if Statmatric*) 


else 
begin 
writeln(SD, ‘Quadrant has less then 2 statistic 
,points *); 
writeln(SD,'and is not used in the overall statics.'); 
end; (*else*) 
end: (*for J*) 


writeln(SD);writeln(SD);writeln(SD); 
writeln(SD, Overall (Real + Imag) statistics for Baud #:',I:4); 


SNROUT: = sqr(BTM[ I] /BTQ[I]) / (BTV[ 1] /(BTQ[ I] -1)); 


SNROUTDB: = 10 * In(SNROUT) / 1n( 10.0); 


writeln(SD,'Overall Baud SNROUT = ':35,SNROUT:10:5,' = ', 
SNROUTDB: 10:5,' DB'); 
writeln(SD,'Overall baud gain (SNROUT / SNRIN) = ':35, 


(SNROUT/SNRIN): 10: 5); 


(*store baud means and variances*) 
SR2| 1):= sk2[ 1] + BIM[ I): 
SR2[2]:= SR2[2] + BTV[I]; 
SR2[ 3]: = SR2[3] + BTQ[I]; 
end; (*for I*) 
writeln(SD);writeln(SD);writeln(SD);writeln(SD); 
writeln(SD, '*** TOTAL OVER ALL ',NUMBAUDS,’ BAUDS ABOVE BAUD ', 
Kx ' : Re yg 
writeln(SD);writeln(SD); 
SNROUT: =sqr(SR2[ 1] /SR2[3]) / (SR2[2] / (SR2[3] - 1)); 
SNROUTDB: =10 * In(SNROUT) / 1n(10.0); 
ueiceln( SD, Gverall SNROUT = '; 35,SNROUT: 10:5,° = i 
SNROUTDB: 10:5,' DB'); 
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writeln('Overall SNROUT = ':35,SNROUT: 10:5,' =', 
SNROUTDB: 10:5,' DB‘); 

writeln(SD,'Overall gain (SNROUT / SNRIN) = ':35, 
(SNROUT/SNRIN): 10: 5); 


writeln(SD);writeln(SD); 
writeln(SD,'Total possible statistical points ',(Kx/8-2)*NUMBAUDS: 3:0, 
Actual number of statistical points ',SR2[3]: 3:0); 


writeln(SD); 

writeln(SD,'BIT ERRORS': 35); 

writeln(SD,'Quadrant 1 = ':30,BITERROR1); 

writeln(SD,'2 = ':30,BITERROR2); 

writeln(SD,'3 = ':30,BITERROR3); 

writeln(SD,'4 = ': 30,BITERROR4); 

BITERROR1: = BITERROR1 + BITERROR2 + BITERROR3 + BITERROR4; 
writeln(SD,'Total Bit errors = ':35, BITERROR1); 
writeln('Total Bit errors = ':30, BITERROR1); 


close(SD); 
end. (“main body*) 
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